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output,  an  inverse  DFT  digital  postprocessor  was  shown  to  yield  the 
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frequency  and  large  time  bandwidth  products  necessary.  An  iterative 
optical  processor  was  then  considered.  This  system  uses  a  linear  LED 
input  array,  2-D  mask  and  1-D  photodetector  array  to  perform  a  vector- 
matrix  multiplication.  With  feedback,  this  system  can  solve  the 
adaptive  radar  equation  using  a  new  iterative  algorithm./  This  system 
is  most  promising,  initial  experiments  have  been  quite  successful,  and 
available  components  make  fabrication  of  such  a  system  attractive. 
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EVALUATION 


Adaptive  processing  for  radar  systems  becomes  complex  as  the  number 
of  adaptive  weights  increase.  Algorithms  used  to  generate  these  weights 
are  conventionally  implemented  through  general  purpose  computers  or  dedi¬ 
cated  hardware  of  which  numerical  operations  are  limited  to  discrete 
parameters  and  the  instruction  rate.  Convergence  time  also  becomes  a 
problem  as  the  number  of  weights  becomes  large. 

In  an  effort  to  assess  the  parallel  channel  potential  of  optical 
processing,  this  effort  was  initiated.  The  number  of  parallel  channels 
approaches  a  continuum  mode  where  optical  focusing  replaces  arithmetic 
operations.  Results  of  the  work  demonstrate  the  feasibility  of  weight 
determination  through  direct  methods  and  closed  loop  techniques. 

Such  approaches  have  the  potential  of  leading  to  efficient  and 


cost  effective  methods  for  realizing  laTge  scale  adaptive  radar  systems. 

VINCENT  C.  VANNICOLA 
Project  Engineer 
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CHAPTER  1 


1.1 


1.  INTRODUCTION 

Adaptive  phased  array  radars  represent  one  of  the  most  demanding  signal 
processing  problems  that  merits  consideration  of  the  use  of  advanced  signal 
processing  techniques  such  as  optical  processing.  In  this  report,  we  summarize 
the  results  of  a  one  year  research  program  in  which  three  different  optical 
techniques  were  considered  for  adaptive  phased  array  radar  processing. 

The  use  of  optical  residue  arithmetic  systems  is  described  in  Chapter  2. 

A  matrix  multiplier  using  residue  arithmetic  is  described  as  are  decimal/residue 
and  residue/decimal  converters.  The  speed  and  parallel  processing  possible  with 
processors  that  operate  in  the  residue  number  system  are  the  main  advantageous 
features  of  such  systems.  An  optical  or  similar  analog  processor  using 
residue  arithmetic  is  likewise  attractive  because  of  the  increased  accuracy 
that  this  type  of  arithmetic  provides.  In  Chapter  2,  we  present  a  new  correla¬ 
tion  approach  that  we  developed  for  the  formulation  of  problems  and  operations 
in  residue  arithmetic.  We  also  present  the  description  of  three  new  optical 
residue  arithmetic  systems  that  perform  decimal/residue  and  residue/decimal 
conversion  and  a  residue  arithmetic  adder.  We  conclude  with  the  design  of  a 
residue  arithmetic  matrix  multiplier  and  our  recommendations  for  future  research 
in  this  area. 

A  modified  coherent  optical  correlator  (COC)  system  for  adaptive  phased 
array  radar  processing  is  described  in  Chapter  3.  In  this  system,  the  inputs 
to  a  multi-channel  correlator  are  the  heterodyne  received  array  signals.  The 
correlation  of  the  received  signals  with  a  reference  array  signal  is  formed  in 
parallel.  When  the  output  of  this  system  is  integrated  along  angular  slits 
using  a  specially  shaped  solid-state  detector,  the  output  is  a  map  of  the 
angular  distribution  of  the  noise  field.  As  we  show,  a  simple  digital  postprocessor 
that  performs  the  inverse  DFT  can  produce  the  desired  adaptive  weights  from  this 


1.2 


output  In  several  restricted  cases.  A  statistical  analysis  of  this  processor 
is  also  provided. 

An  iterative  optical  processor  (IOP)  system  using  a  linear  input  LED 
array,  2-D  mask  and  1-D  output  detector  with  feedback  is  then  described  in 
Chapter  4.  This  system  performs  a  vector-matrix  multiplication  and  with  feed¬ 
back  results  in  an  iterative  processor  that  has  the  covariance  matrix  as  an 
input  and  that  solves  the  vector-matrix  adaptive  phased  array  radar  equation 
for  the  desired  weights  using  a  new  iterative  algorithm. 

For  all  three  candidate  optical  processors,  laboratory  optical  systems 
were  assembled,  tests  and  demonstrations  of  each  were  conducted,  and  analysis 
was  conducted  of  the  system's  key  features,  a  study  of  the  availability  of 
components  was  made  and  recommendations  for  future  work  are  advanced. 

1.2  ADAPTIVE  PHASED  ARRAY  RADAR  PROCESSING  REVIEW 

A  simplified  block  diagram  of  an  adaptive  phased  array  processor  is  shown 
in  Fig.  1.1.  The  basic  concept  is  to  multiply  each  of  the  received  signals  v 
by  an  appropriate  weight  and  to  then  sum  these  products  to  produce  the  output 

N 

E  -  l  Wn  V  (l.D 

n 

If  the  same  weight  is  applied  to  all  elements,  the  beam  formed  is  normal  to  the 
•tray  and  described  by 


E  (fl) 


,  ,tt  Nd  r.ln  0, 
Bin  ( - r - ) 


,7i  d  sin  0. 
Sin  ( - r - ) 


(1.2) 


where  d  is  the  element  element  spacing,  0  is  the  steering  vector  of  the  array, 
•nd  N  is  the  number  of  elements.  Dy  varying  the  phases  of  the  weights  w  ,  the 


beam  can  be  steered  to  different  direct'ons.  This  is  achieved  by  selecting  the 
phases  such  that  a  signal  incident  from  the  desired  scan  angle  adds  in  phase 
across  the  array. 


By  adjusting  the  amplitudes  and  phases  of  the  w^,  the  side  lobe  levels  can 

be  decreased  and  the  effects  of  other  noise  sources  in  the  antenna's  field  of 

view  can  be  reduced.  In  such  APAR  systems,  nulls  are  placed  at  angles  and 

frequencies  in  the  antenna  pattern  corresponding  to  different  noise  sources. 

This  adaptive  control  is  achieved  by  separate  adaptive  loops  on  each  of  the 

antenna  elements  as  shown  in  Fig.  1.2.  The  steering  signals  s*  indicate  the 

direction  in  which  the  array  is  steered  and  hence  the  location  of  the  antenna's 

main  beam.  The  ur  values  are  the  correlation  of  the  received  array 

•signals  v  and  the  output  g  «*  Y  w  v  .  When  w  (*)  g  is  large,  u  and  hence  w 
b  n  L  n  n  n  n  n 

change  rapidly.  The  purpose  of  an  adaptive  array  is  to  reduce  the  noise  in  g. 

When  u^  in  a  given  channel  is  large,  the  corresponding  weight  w^  has  a  larger 

effect  on  reducing  the  residue  noise  in  g. 

For  receiver  noise  only,  the  same  for  all  channels,  the  weights  v  all 

approach  the  same  value.  If  one  vn  is  larger,  the  corresponding  w^  v^  term 

in  g  is  larger  and  this  will  cause  larger  and  u^  values  for  that  channel, 

which  will  decrease  w^.  Thus,  an  adaptive  array  achieves  uniform  illumination 

only  for  receiver  noise  alone.  When  the  received  signal  energy  is  less  than 

the  interference  and  noise  energy,  the  adaptive  loops  attempt  to  minimize  the 

input  power  (subject  to  the  steering  vector  constraint).  If  the  noise  in  g 

is  spproximately  zero,  the  z  values  will  be  small  and  the  w  will  be  rcla- 

n  n 

tlvely  constant. 

The  weights  can  be  described  by 


G 
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The  low  pass  filter  with  time  constant  t  can  be  described  by 


r  u  +  u 
n  n 


(1.4) 


and 


z 

n 


V*  g 


(1.5) 


These  equations  describe  the  responseof  each  of  the  N  adaptive  loops.  The 
variation  of  w^  depends  on  the  inputs  from  all  N  channels  and  on  the  weights 
of  the  other  N  -  1  channels.  In  matrix  form,  we  describe  the  loops  by 


W  -  G  (S*  -  U) 


T  U  +  U 


Z 


Z  «  V  *  (V  -  W) , 


(1.6) 


where  V^,  is  a  row  vector  and  W  etc.  are  column  vectors.  In  terms  of  average 
values 


*  T 

Z  -  V  V  W  -  MW, 


(1.7) 


where  M  is  the  covariant  matrix 

Mjj  »  Avg  (V*  Vj). 


(1.8) 


If  ft  la  a  constant  steering  vector,  U  -  W/C  and  the  loop  is  described  by.  . . 


T  V/G  +  (M  +  I/G)  W  -  S*.  (1.9) 

vbaro  I  is  the  identity  matrix.  Equation  (1.9)  describes  a  set  of  N  different 
equations  that  describe  the  average  response  of  the  array  weights  W.  The  values 
kl  depend  upon  the  external  noise  (thru  M)  and  the  control  loop  parameters  (S*, 

C  sod  t).  The  conventional  approach  is  to  solve  (1.9)  for  W  and  then  apply  these 
to  the  received  array  signals. 

The  covariance  matrix  M  is  fundamental  to  all  APAR  processing  theory  since 
It  describes  the  noise  environment.  The  diagonal  terms  in  M  are  a  measure  of 
the  power  in  each  channel,  whereas  the  off  diagonal  elements  of  M  describe  the 
dir  ction  of  arrival  of  the  noise.  The  rate  of  convergence  of  (1.9)  depends 
upon  the  noise  environment.  The  steady  state  solution  to  (1.9)  is  found 
with  W  ■  0  (assuming  G  >>  1)  to  be 


Wrs*-  (1.10) 

The  desired  weights  W  can  then  be  found  from  S  and  M  *  as 

W  -  M_1  S*.  (1.11) 

Equation  (1.11)  can  be  obtained  by  inverting  the  matrix  M,  and  forming  the 
indicated  vector/r.atrix  product.  Computation  of  M  requires  2  N  samples  at 
the  signal  frequency  to  yield  acceptable  statistics. 

Various  correlation  loop  processors  exist.  The  maximum  SNR  circuit  is  the 
most  popular.  It  uses  the  Widrow  least  mean  square  circuit.  It  requires  an 


Initial  estimate  for  W,  computation  of  the  gradient  at  different  surface  points, 
etc.  until  the  minimum  concave  surface  is  obtained.  The  system  In  Fig.  1.2 
uses  a  typical  Hovells-Appelbaum  loop  in  which  the  residue  is  fed  back  and 
correlated  with  each  of  the  received  signals.  The  filtered  outputs  are  pro¬ 
portional  to  the  gradient.  When  subtracted  from  the  steering  signal,  they  yield 
the  adaptive  weights.  A  modified  random  search  technique  with  many  variations 


CHAPTER  2 


A  RESIDUE  ARITHMETIC  PROCESSOR  FOR 
ADAPTIVE  PHASED  ARRAY  RADAR 


2.1  INTRODUCTION 

The  data  processing  requirements  for  adaptive  phased-array  radar  (APAR) 
systems  represent  one  of  the  largest  challenges  to  present  technology.  For 
this  reason,  new  and  fresh  approaches  to  APAR  processing  are  necessary 
to  handle  the  massive,  data  processing  requirements  foreseen  for  the  future. 

The  novel  data  processing  approach  considered  in  this  chapter  is  the  use  of 
residue  arithmetic  (RA) . 

Addition,  subtraction  and  multiplication  can  be  performed  in  RA  with  no 
carries  required.  Because  of  this,  RA  is  directly  implementable  in  a  parallel 
processor  and  is  thus  one  of  the  fastest  computational  methods  available. 

Since  an  optical  system  is  inherently  a  parallel  processor,  a  RA  optical  pro¬ 
cessor  implementation  is  most  attractive.  RA  is  also  attractive  because  it 
offers  a  solution  to  the  problems  of  limited  accuracy  and  limited  dynamic  range 
present  in  many  optical  and  analog  data  processors. 

A  brief  review  of  residue  arithmetic  is  included  in  Sect.  2.2.  The  rea¬ 
sons  for  considering  this  novel  processing  approach  (speed,  parallel  processing 
accuracy,  high  dynamic  range,  etc.)  can  then  be  seen  more  clearly  (Sect.  2.3). 
The  first  major  issue  to  be  decided  upon  in  a  residue  processor  is  how  to 
realize  the  cyclic  permutations  required  to  represent  numbers  in  moduli.  In 
Sect.  2.4,  we  discuss  the  various  alternatives  available  and  the  reasons  why 
we  chose  pulse-position  coding  to  represent  numbers  in  RA. 

The  second  issue  to  be  addressed  in  a  RA  processor  is  the  formulation  of 
the  RA  operations  performed.  One  of  the  major  original  achievements  of  this 
contract  was  the  correlation  formulation  of  RA  operations.  This  is  included 
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in  Sect.  2.7.  Specific  RA  system  designs  for  decimal-to-residue  and  residue- 
to-decimal  converters  are  provided  in  Sects.  2.5  and  2.6.  As  part  of  this 
program,  we  provided  the  first  experimental  demonstration  of  an  optical  RA 
system  (Sect.  2.9).  The  design  of  a  RA  adder  (Sect.  2.8)  shows  how  the  pro¬ 
cessing  system's  architecture  directly  follows  from  our  correlation  formulation 
approach.  Another  advantage  of  our  correlation  formulation  of  RA  is  the  ease 
with  which  RA  processors  can  be  realized  using  many  diverse  technologies 
(LSI,  CCDs,  SAW,  Acousto-Optics,  Waveguides,  Integrated  Optics,  etc.). 

In  Sect.  2.10,  we  quantify  the  dynamic  range  obtainable  in  a  RA  optical 
processor  and  suggest  several  added  schemes  to  increase  the  system's  dynamic 
range.  A  related  aspect  of  a  RA  optical  processor  that  emerges  clearly  and 
numerically  in  Sect.  2.10  is  the  utilization  of  the  vast  space  bandwidth 
product  available  in  an  optical  processor. 

In  Sect.  2.11,  we  briefly  review  another  approach  to  the  realization  of 
residue  operations  (the  use  of  maps).  We  also  include  system  realizations  of 
residue  processors  using  technologies  such  as  ROMs  and  integrated  optics.  The 
design  of  a  multiplier  in  residue  is  then  discussed  in  Sect.  2.12.  With  these 
basic  building  blocks  (converters,  adders  and  multipliers),  we  sketch  the 
design  of  a  parallel,  pipelined  matrix  multiplier  using  residue  arithmetic 
in  Sect.  2.13.  Conclusions  are  then  advanced  in  Sect.  2.14. 

2.2  RESIDUE  ARITHMETIC 

A  brief  review  of  residue  arithmetic  follows.  In  the  residue  number 
system,  an  integer  is  represented  by  the  N-tuple  set  of  residues 
(R,,,  .Rn,^,  ...  Rmjj)  with  respect  to  the  N  different  relatively  prime  integer 
moduli  m^,  ...  ,  m^.  The  residue  Rm  is  the  least  positive  integer  remainder 
of  the  quotient  of  X  and  ra^.  The  maximum  integer  value  that  can  be  represented 
by  N  moduli  is  M-l  where 
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i 


* 


N 

M  =  H  m  . 
i  =  l  1 

An  example  will  easily  clarify  this.  To  represent  13  in  the  residue 
system  with  moduli  (5, 7, 9, 2),  we  divide  13  hy  the  first  modulus  5  and  obtain 
2  and  a  remainder  of  3.  Dividing  13  by  the  second  modulus  7,  we  obtain  1 
and  a  remainder  of  6.  Dividing  13  by  =  9,  the  remainder  is  4.  For  m^  =  2, 
the  remainder  is  1.  These  four  remainders  are  the  representation  of  13  in 
the  residue  number  system  with  moduli  (5, 7, 9, 2),  i.e.  13  =  (3, 6, 4,1). 

To  see  the  advantages  of  residue  arithmetic,  we  consider  how  to  realize 
the  sum  of  13  and  59  in  residue  with  moduli  (5, 7, 9, 2).  Proceeding  as  before, 
we  find  59  *  (4, 3, 5,1).  To  sum  two  numbers  in  residue  notation,  we  sum  each 
separate  pair  of  residue  numbers  in  the  same  modulus  iik  and  retain  only  the 
residue,  i.e. 

m.  =  (5, 7, 9, 2) 

11  =-(3,6,4,17 
+59  =  (4,3,5,!) 

72  =  (2, 2, 0,0) 

where  3+4  modulo  5  has  a  remainder  2,  1*1  modulo  2  is  0,  etc.  Checking,  we 
find  that  the  decimal  sum  72  is  (2, 2, 0,0)  in  moduli  (5, 7, 9, 2). 

To  subtract  two  numbers  in  residue,  e.g.  72-13,  we  invert  13  in  residue 
(bv  simply  complementing  each  bit  i  modulo  nu  in  the  residue  representation) 
i.e. 

m.  =  (5, 7, 9, 2) 

13  =-(776,4,1) 

-13  =  (2, 1,5,1) 

The  bits  i  of  the  residue  numbers  72  and  -13  are  then  added  as  before 

72  =  (2, 2, 0,0) 

-13  =  (2, 1,5,1) 

59  =  (4, 3,5,1) 
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To  multiply  two  numbers  in  residue,  we  convert  each  number  to  its  residue 
representation,  multiply  each  bit  pair,  and  retain  the  residue  of  each  product 
modulo  m^.  Consider  forming  the  product  of  19  and  12  in  the  same  moduli. 

">i  =  (5, 7, 9, 2) 

19  =  (4, 5, 1,1) 

12  -  (2, 5. 3.0) 

19  x  12  =  (3, 4, 3,0)  =  228 

where  1  times  0  is  0  in  any  modulo,  3  times  1  is  3  modulo  9,  5  times  5  is 
25  but  modulo  7  the  remainder  or  residue  is  4,  and  4  times  2  is  8  which  yields 
a  residue  of  3  in  modulo  5.  The  product  of  19  and  12  in  decimal  is  228.  But 
228/5  equals  45  and  a  residue  of  3,  etc.  from  which  228  =  (3, 4, 3,0)  with 
respective  moduli  (5, 7, 9, 2).  From  the  abova  examples,  the  attractive  features 
of  residue  arithmetic  are  apparent. 

2.3  SPECIAL  FEATURES 

As  shown  in  Sect.  2.2,  no  carries  are  required  in  a  residue  arithmetic 
computation.  Thus  this  is  inherently  one  of  the  fastest  possible  computing 
methods  available  [2.1  -  2.3]  and  thus  of  potential  concern  to  the  matrix 
inversion  process.  A  second  feature  of  residue  arithmetic  is  that  in  residue 
processing  a  given  computation  is  divided  into  subcomputations  of  reduced 
complexity.  This  allows  these  sub-operations  to  be  performed  in  parallel 
with  reduced  dynamic  range  requirements.  The  accuracy  of  the  resultant  opera¬ 
tions  is  high  because  the  system's  dynamic  range  is  proportional  to  the  product 
of  the  residue  moduli  used.  Thus,  residue  arithmetic  directly  allows  the  use 
of  multi-channel  parallel  computations  that  are  the  forte  of  optical  computing 
and  overcomes  any  potential  dynamic  range  and  accuracy  limitations  often 
associated  with  optical  processors. 
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2.4  REALIZATIONS 

Residue  processing  requires  a  cyclic  permutation  procedure.  This  can  be 
realized  by:  polarization,  phase,  or  pulse  position  coding.  The  former  two 
approaches  are  being  pursued  by  Collins  [2.4].  It  appears  unrealistic  in 
practice  since  nearly  ideal  2-D  spatial  light  modulators  are  required.  Since 
such  devices  are  space  variant,  we  considered  pulse  position  modulation. 

Huang  [2.5  -  2.6]  was  the  first  to  suggest  an  optical  residue  arithmetic 
processor.  Recent  work  [2.7  -  2.8]  has  concentrated  on  similar  pulse  position 
coding.  We  concur  with  such  an  approach  because  of  the  large  overall  error  intro¬ 
duced  into  a  residue  processor  when  an  error  occurs  in  any  bit. 

2.5  DECIMAL-TO-RESIDUE  CONVERSION  [2.9] 

We  consider  the  conversion  of  a  decimal  number  into  the  residue  number 
Rj,^  modulo  and  the  specific  correlator  topology  in  Fig.  2.1.  We  define  a 
unit  decimal  distance  Ax  and  thus  represent  the  decimal  number  x  by  a  delta 
function  at  xQ  =  XAx  in  Pq.  The  amplitude  transmittance  of  Pq  is  thus 


g(x0)  =  iS(x0-XAx). 


(2.1) 


Lenses  L^  in  Fig.  1  form  the  1-D  Fourier  transform  (FT)  of  g  horizontally , 
while  imaging  vertically  (i.e.  a  multi-channel  1-D  FT).  We  consider  only 
one  y  channel  for  simplicity.  The  FT  of  g  incident  on  is 

C(u)  =  exp  (j2n  uXAx) .  (2.2) 

Spatial  distance  x^  in  1’  is  related  to  input  spatial  frequency  u  by 


X1  =  XfLU’ 


(2.3) 


r 
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where  \  is  the  wavelength  of  the  laser  light  used  and  fj  =  f  ^  is  the 
focal  length  of  lens  in  Fig.  2.2  Eq.  2. 2)describes  a  planewave  incident 
on  at  an  angle  that  depends  on  XAx  (the  position  of  the  input  delta  function 
and  hence  the  value  of  the  decimal  input  number  X) . 


To  realize  the  conversion  of  X  to  modulo  we  place  at  P^  a  square 
wave  grating  with  transmittance 


11.  (u)  =  y  exp  (j2nunm.Ax),  (2.4) 
n  = 

where  the  fundamental  frequency  ut.  of  the  grating  present  on  this  channel  i 
is  chosen  such  that  the  separation  between  the  dc  and  +  first  order  terms  in 
the  FT  of  (  2.4)  is 


m.  Ax  =  u  .if,, 

r  i'.t  1. 


(2.5) 


For  simplicity,  the  same  unit  distance  Ax  is  assumed  for  decimal  nails  in  1’^ 
and  for  residue  units  in  P^  (this  is  equivalent  to  assuming  a  1:1  imaging 
system  for  and  l.^ ) . 

In  one  version  of  a  decimal/residue  converter,  the  input  Pq  data  is  a 
vertical  slit  (one  decimal  number)  and  the  desired  outputs  at  P2  are  the 
residue  numbers  R_  for  N  moduli  m.  as  shown  in  Fig.  2.1.  In  such  a  case, 

ui^  X 

the  system  will  have  N  channels  and  at  P^  there  will  be  N  different  square 
waves  recorded  on  these  N  channels.  The  frequency  u  ^  of  the  grating  on 
channel  i  will  correspond  to  the  modulus  m^  of  that  channel  as  in  (2.5). 
Again,  for  simplicity,  we  restrict  attention  at  present  to  one  channel  i. 

The  light  amplitude  distribution  leaving  channel  i  of  plane  1’^  is  the 
product  of  (2.2)  and  (2.4)  or 

HU .  -  }  exp  [  j2iiu(X-nm.)Ax] .  (2.6) 

n  = 


Lenses  L ^  in  Fig.  2.1  form  the  1-D  horizontal  FT  of  (2.6)  and  af  we  find 

t2(-K2)  =  £  >M\2-(X-nm.  )Ax]  ,  (2.7) 

n  = 

where  x?  is  the  spatial  coordinate  of  P,.  Tliis  pattern  consists  of  delta 
functions  with  center-to-center  spacings  in^Ax  proportional  to  the  modulus  nu 
and  replicated  at  distances  that  are  multiples  of  this  modulus.  These  delta 
function  outputs  occur  at 


x2  =  (X-min) Ax  .  (2.8) 

The  pulse  position  coding  of  X  in  1’^  .and  Lhe  carrier  modulation  [multipli¬ 
cation  of  (2.2)  by  (2.4)  and  Fourier  transforming  the  result]  operations  involved 
have  now  been  explained.  Before  describing  the  last  required  operation 
(aperture  control),  we  digress  for  the  moment  to  a  general  description  of 
decimal/residue  conversion.  Any  decimal  number  X  can  be  written  as 


where  n  is 
or  residue 


X  =  nm.  +  R  , 

l  m. 

l 


the  number  of  times  X  is  divisible  by  nn 


R  must  satisfy 
m . 

l 


(2,9) 

and  where  the  remainder 


0  <  K  in .  -  1  . 

—  i 

i 


(2.10a) 


Solving  (2.9)  for  R  ,  we 


l;rom(  2. 11a)  ,  we  see  that  R 


f  ind 


R 

m . 

t 

is  the 


X-nm .  . 
i 

remainder 


after  m. 

l 


(2.11a) 


has  been  subtracted  n 


times  from  X.  Multiplying  (2.10a)  and  (2.11a)  by  Ax,  we  obtain 
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Fig.  2.1 


Schematic  diagram  of  a  decimal-to-residue  optical  converter 
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’ig.  2.2  Schematic  diagram  of  a  residue-to-decimal  optical  converter, 
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0  <  R  Ax  <  (m.-l)Ax  (2.10b) 

—  m.  -  l 
l 

R  Ax  =  (X-nm.)Ax.  (2.11b) 

in .  ) 

l 

Now,  returning  to(2.7)Jwe  see  tliat  tlie  locations  of  the  n  output  delta 
functions  in  P ^  of  Pig.  2.1  satisfy  (2.8),  which  agrees  with  (2.11b).  To  insure 
that  the  subtraction  of  m  from  X  has  been  performed  the  proper  number  of 
times  n  so  that(2.10a)  and  (2.10b)  are  satisfied,  we  examine  only  the  P„  region 

satisfying 

0  <_  x.  <  (ra.-i)Ax.  (2.12) 


Tills  is  obviously  equivalent  to  satisfy ing(2 . 10a)  and  (2. 10b). Only  one  of  the 
delta  function  outputs  in  (2.8)  satisfies  (2.12). 

If  we  place  a  rectangular  aperture  at  of  Fig.  2.1  of  width  m^Ax  in  x2 
centered  at  P9  as  defined  by (2. 12) , t here  will  be  tie  possible  locations  of  the  out¬ 
put  spot  of  light  witlri  ibis  aperture.  The  location  of  the  output  peak  of  light 

within  this  aperture  denotes  the  residue  K  modulo  m.  of  the  decimal  input 
1  m.  l 

i 

X  in  terms  of  a  unit  residue  distance  Ax  (i.e.  pulse-position  coded).  By 
aperture  control  we  simply  mean  placing  an  aperture  at  P.,  to  select  the 

remainder  or  residue  in  the  proper  range. 

The  above  analysis  considered  only  one  channel  of  the  system  of  Fig.  2.1. 

If  this  system  is  to  be  used  to  convert  the  decimal  number  X  into  the  N 

residues  1’^  for  the  N  moduli  ng  ,  the  input  I’  pattern  is  a  slit  as  shown  in 
i 

Fig.  2.1. Its  Xy  location  XAx  denotes  the  decimal  input  number  X.  At  P^,  we 

now  place  N  gratings  of  spatial  frequencies  u^ .  (corresponding  to  the  N 

moduli  tin  to  be  used  to  represent  X)  on  N  separate  lines  or  channels  at  P^. 

Within  Lhe  aperture  placed  at  1’2  there  will  be  N  channels  with  an  output  spot 

of  light  at  a  different  horizontal  location  on  each  channel.  The  horizontal 

location  ol  the  output  on  each  channel  i  denotes  the  residue  number  R  for  X 

in. 

l 
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in  the  modulus  m±  for  the  particular  channel.  This  particular  arrangement 
of  a  decimal-to-residue  converter  may  not  be  optimum  from  a  system  standpoint. 

We  could  also  insert  various  decimal  numbers  on  separate  channels  at  PQ  in 
Fig.  2.1.  By  use  of  a  single  grating  at  P2,  the  conversion  of  N  decimal 
inputs  into  their  equivalent  residue  representations  with  respect  to  one 
modulus  results  at 

2.6  RESIDUE/DECIMAL  CONVERSION 

The  schematic  of  a  residue-to-decimal  optical  converter  is  shown  in  Fig.  2.2. 

The  system  is  topologically  the  same  as  that  of  Fig.  2.1,  however  the  inputs  at  PQ 

are  now  N  residue  numbers  in  the  N  moduli  m.  chosen.  These  N  input  residue 

numbers  R  arc  represented  by  a  pulse  of  light  whose  horizontal  distance  RAx 
i 

along  channel  n  represents  the  residue  number  by  pulse-position  coding  as 
before. 

At  1'^  of  fig. 2.2, N  gratings  are  placed  on  the  N  diflereul  channels  with 
the  spat ia l  frequency  of  each  grating  chosen  to  correspond  to  Lbe  modulus 
tin  for  that  channel.  The  output  along,  a  given  channel  n  at  I’  of  Fig.  2.2 
will  be  a  set  of  della  functions  as  before.  The  horizontal  posit  ions  ol 
those  delta  functions  now  correspond:;  Lo  all  possible  decimal  numbers  to 
which  the  input  residue  number  could  correspond  (for  the  particular  modulus 
of  that  channel).  The  horizontal,  position  XAx  in  I’7  at  which  the  delta  func¬ 
tions  on  all  N  rhanne  1  s  are  i  l  i  p.ned  *•*:  Hu-  proju*r  ilrrini.il  mu  pul  \  ror- 

respondinj;  to  the  N  input  residue  numbers  K  in  the  N  moduli  m .  ,  To  detect 

m  .  i 

i 

the  output.,  a  set  of  hori  /.onLa  l  1  y-oriented  linear  photo  diode  or  c:r,l)  self- 
scanned  detector  arrays  can  he  placed  at  P  .  The  output  from  an  N  input 
NOR  or  AND  pule  (depend  i  up,  on  the  polarity  outputs  desired)  in  svinh  run  i.  y.,\ t  i  on 


m- 
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with  a  Ax  couiiLor  will  denote  the  decimal  output  value  X  (or  equivalently 
the  horizontal  X.\x  position  in  aL  which  all  delta  functions  on  a L J 
channels  .ne  present  1 . 

The  optical  realization  suggested  in  Fig.  2.2  is  similar  in  principle 

to  a  scheme  suggested  earlier  by  Huang  [2.6],  In  the  version 

in  Ref.  2.6  (Fig.  6),  input  light  was  passed  sequentially  through  gratings 

(corresponding  to  each  m^) .  However  in  this  prior  scheme  the  lateral  posi- 


achieve  the  proper  superposition  at  the  output.  In  our  scheme  suggested  in 
Fig.  2.2  mechanical  motion  of  the  gratings  is  not  required. 

2 .7  CORRELATION  FORMULATION 

The  following  general  correlation  formulation  of  residue  operations  will 
be  more  clear  now  that  detailed  descriptions  of  the  two  converters  have  been 
provided.  The  designs  in  Figs.  2.1  and  2.2  will  also  be  clearer  from  the 
following  formulation.  In  fact,  an  attractive  feature  of  this  correlation- 
based  description  is  that  the  system  architecture  follows  directly  from  the 
analytical  formulation.  Such  a  correlation  formulation  is  most  general, 
since  it  allows  realization  of  residue  operations  using  other  technologies. 

We  have  selected  a  coherent  optical  realization  because  it  is  one  of  the 
easiest  to  demonstrate  and  because  of  the  real-time  and  parallel-processing 
advantages  of  such  processors. 

We  consider  the  dec imal-to-residue  conversion  of  X  to  R  modulo  m..  Now 

m .  r 

r 

that  the  pulse-position  coding,  carrier  modulation  and  aperture  control 

concepts  have  been  described  ,  the  general  description  of  the 

system  as  a  correlator  is  straightforward.  The  input  is  a  decimal  number  X 

represented  by  g(xQ)  =  <$(x()-XAx).  The  desired  residue  output  R^  is 

f(x^)  =  6  [x^- (X-nnu  )  Ax  ]  .  The  required  system  impulse  response-  needed  to 


convert  g  into  f  is 
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U(Xq)  =  £  5 (xQ-nnuAx) ,  (2.13) 

n 

where  h  is  simply  the  FT  of  the  grating  at  described  by  (2.4). 

Since  the  system's  impulse  response  is  real  and  symmetric,  the  convolution  and 
correlation  operations  arc  equivalent  and  the  system's  output  at  1’2  is  the  correla¬ 
tion  of  g  and  h 

g  ©  h  =  )]  ft(x()-X Ax)  iS (Xy-nuuAx  +  x^)dx0 

n 

=  l  /6(xq+x2)  <S  [ Xq-  (X-iutk  )  Ax  ]  dxQ 
n 

=  J  6[x2  -  (X-nm.)Ax),  (2.14) 

n 

The  minimum  positive  location  x2  of  the  output  delta  function  in  the  ?2 

interval  0  <_  x2  _<  (m^-l)Ax  (where  x2  is  now  the  shift  parameter  in  the 

correlation  or  equivalently  the  output  coordinate)  is  the  desired  residue 

output  .  Residue-to-decimal  conversion  can  be  formulated  similarly. 

The  general  system  block  diagram  for  a  residue  converter  is  shown  in 

Fig.  2.3.  The  input  is  pulse-position  coded  by  a  deflector  or  other 

component,  the  correlator  can  be  realized  using  many  technologies,  the 

output  window  function  is  equivalent  to  looking  only  in  the  proper  region 

of  output  correlation  space.  The  approach  allows  one  to  more  easily 

realize  a  residue  arithmetic  system  in  terms  of  the  many  technologies  that 

2 

use  correlators  (CCD,  SAW,  I  L,  Integrated  optics,  etc.).  We  consider 
only  coherent  optical  techniques  in  the  specific  system  designs  to  be 
described,  thereby  utilizing  the  space  bandwidth,  parallel  processing  and 
real  time  throughput  of  such  systems.  However,  the  correlation  formulation 
presented  is  general  enough  to  allow  residue  systems  to  be  realized  using 
other  technologies.  The  specific  optical  system  architectures  presented 
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Fig.  2.3  General  correlator  schematic  for  a  decimal /residue  converter. 


Fig.  2 .4  Schematic  diagram  of  a  residue  arithmetic  optical  adder. 
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and  experimentally  demonstrated  have  been  included  to  better  convey  the 
operations  required  and  as  examples  of  transitions  from  theory  to  system 
architecture  to  experimental  verifications. 


2.8  RESIDUE  ARITHMETIC  OPTICAL  APDUIl 

To  demonstrate  the  utility  of  the  correlation  formulation  of  residue 

arithmetic,  we  describe  the  addition  of  the  two  residue  numbers  R  and  R. 

a  b 

modulo  nu  as  a  correlation.  We  then  sketch  the  design  of  a  residue  arithmetic 
adder  based  on  this  formulation. 

The  two  numbers  to  be  added  are  represented  in  pulse-position  coding  by 


a  =  6 (x  -R  Ax)  (2.15a) 

()  a 

b  =  6(x  -K,  Ax)  (2.15b) 

U  b 


The  desired  sum  is 


where  the  subscript  itk 
modulo  nu.  To  produce 
the  mirror  image  of  b, 

a  0  b  =  / 

=  f 


Ik  +  K.  I  (2.16) 

'a  b  ui . 

i 

denotes  that  the  output  ol  the  sum  is  the  residue 
the  desired  result,  we  torm  the  correlation  ol  a  and 

<S  (x-R^Ax)  6  (x+R^Ax+x*  )dx 

£ 1  x-  (R  +R.  )  Ax  ]  6  (x  1-x '  )dx 
a  b 


.Slx’-CK  +Rb)Ax) 


(2.17) 


Prom  this,  we  see  that  the  output  of  the  above  correlation  is  a  delta 
function  located  in  correlation  space  at  x'  =  (R  +Rb)Ax  where  x’  is  the  correla¬ 
tion  shift  parameter.  This  corresponds  to  the  algebraic  sum  of  the  two  residue 
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numbers  pulse-posit ion  coded.  To  convert  this  sum  to  the  residue  modulo  m^,  we 

use  a  grating  and  output  aperture  as  in  Fig.  2.1  to  produce  the  desired  result  in 
(2.16). 

The  schematic  of  an  optical  system  to  realize  these  operations  is  shown  in 
Fig.  2.4. The  two  pulse-position  coded  input  residue  numbers  to  be  added  arc 
placed  side-by-side  in  the  input  plane  I’q  of  a  joint  transform  correlator.  We 
describe  the  transmittance  of  P^  in  more  conventional  optical  notation  (in  1-D 
for  the  case  of  one  modulus  only  for  simplicity)  by 

W  =  4(x0"xa)+6(x0+xb)’  (2.18) 

where  x  =  R  Ax  and  x.  =  R,  Ax  denote  the  locations  of  the  two  delta  functions 
a  a  b  b 

(for  R^  and  R^)  in  P^.  Lens  forms  the  1-D  Fourier  transform  at  of  t^ 

in  the  x  direction  while  imaging  in  the  y  direction. 

The  light  distribution  in  1-D  incident  on  P^  is  thus 

f^u)  =  cxp(j2nuxa)+exp(-j2iruxb)  .  (2.19) 

An  optically  addressed  spatial  light  modulator  [2.10]  (SLM)  is  placed 

2 

at  P1>  The  transmittance  of  P^  after  exposure  is  | f ^ |  or 

j2nuxa  -j2nuxb  ^ 
t^  =  e  +  e 

=  2  +  2  cos  [2nu(x. +x,  )  ] .  (2.20) 

3  D 

This  SLM  is  read  in  reflection  from  the  beam  splitter  BS  in  Fig. 2.4.  Lens  L£ 
in  Fig.  2.4  images  P^  onto  At  >  we  Place  a  grating  whose  transmittance  in 

1-D  for  the  y  channel  of  concern  we  represent  by 
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-j  2trunm .  Ax 

c2  =  l  e  (2.21) 

n 

Tftfe  light  distribution  leaving  is  now 

=  l  exp  [j2nu(x  +x. -nm.Ax) J.  (2.22) 

n  3  t>  i 

Lens  L3  iw  FiS-  2.4  forms  the  FT  of  (2.22)  at  P-j.  This  P3  pattern  is  a  set  of 
delta  functions, 

f3  =  l  5tx3  -  (xa+xh-nm.Ax) ]  (2.23) 

n 

=  I  ^  t x n  -  (U  +K.  -um  .  )  Ax  ]  . 

O  ci  1)  1 

n 

The  location  of  the  delta  function  at  P.^  of  Pig.  2.4  in  the  range 

0  1  x3  1  (mi-l)Ax  provides  the  desired  residue  of  the  sum  of  the  two  residue 

numbers  moduio  m.. 

l 

In  practice  a  and  l>  are  represented  bv  N  residue  numbers  (Cor  the  case 

of  N  moduli  nu).  One  way  in  which  the  contents  of  I’y  in'  l'ig.  2.4  can  be 

arranged  to  add  two  residue  numbers  with  N  moduli  is  to  allot  N  channels  in 

each  half  of  plane  P_.  The  N  pairs  of  residue  numbers  R  and  R,  for  the 

0  am .  bm . 

1  i 

N  moduli  can  be  recorded  on  these  N  separate  channels  in  pulse-position  code. 

For  this  input  data  arrangement,  we  would  use  N  gratings  at  N  dilioront 
spatial  frequencies  (corresponding  to  the  N  moduli)  on  N  lines  at  I’,,  of  l  ig.  2.4. 

At  of  Fig. 2/4, there  will  also  lie  N  channels  with  the  position  of  the  delta 
function  on  each  vertical  channel  proportional,  to  the  residue  of  the  sum  of 
the  corresponding  residue  numbers  in  the  corresponding  modulus  of  that  channel. 

The  adder  element  shown  in  Fig. 2. 4  is  one  of  the  basic  elemental  building  blocks 

of  a  residue  arithmetic  system  (a  multiplier  can  be  realized  by  successive  additions^ 
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Tlie  system  of  Fig.2.1was  assembled  to  demonstrate  and  verify  the  decimal/ 
residue  conversion  method  described  earlier. 

To  keep  the  unit  decimal  distance  Axfi  in  and  the  unit  residue  distance  Ax£ 
in  1*2  equal  to  Ax,  a  1:1  imaging  system  was  used.  The  focal  lengths  of  both 
cylindrical  lenses  were  f  =  300  mm  witli  f |  =  762  mm  focal  lengths  for  both 
spherical  lenses. 

As  noted  several  times  earlier,  the  optimum  decimal/residue  converter 
need  not  always  convert  a  decimal  input  number  into  its  N  residue  numbers.  To 
demonstrate  an  alternate  scenario  and  to  still  demonstrate  the  basic  principal 
of  decimal/residue  conversion  as  a  correlation,  the  21  decimal  numbers  0  to  20 
were  used  as  the  inputs  on  21  channels  in  of  Fig. 2.1.  We  consider  the  con¬ 
version  of  these  21  decimal  inputs  to  21  residue  numbers  all  in  the  same  modulus. 
The  input  used(Fig.  2. 5a) consists  of  an  impulse  on  each  of  the  21  vertical 
channels  in  1’q.  The  horizontal  position  of  each  of  these  impulses  denotes  a 
different  decimal  input  number.  Since  the  sp.icings  between  impulses  on  suc¬ 
cessive  channels  are  equal  to  Ax  =  0.7  mm  (the  unit  input  decimal  distance), 
the  decimal  input  on  the  bottom  channel  0  Is  0,  the  number  on  channel  1  is  1, 
etc.,  up  to  the  top  channel  which  represents  the  decimal  number  20  in  pulsc- 
posiliou  code. 

This  input  pattern  was  produced  by  superimposing  a  titled  slit  over  a  1.4 
cy/mm  grating.  The  input  pulses  are  then  separated  by  1/1.4  =  0.7  mm,  =  Ax. 

Since  the  duty  cycle  of  the  1.4  cy/mm  grating  was  only  6%,  the  width  of 
the  input  pulse  is  only  0.06  (0.7)  mm  =  40  pm  or  far  less  than  the  0.7  nm 
decimal  limit  used.  The  capacity  of  the  system  is  thus  not  even  approached 
by  the  example  shown. 


At 


of  Fig.  1,  we  place  a  grating  of  spatial  frequency  u  =10  cy/mm. 


From(2.5),we  see  that  Ax  and  must  satisfy 


Ax  =  0.7  mm  =  u  Xf ,/m. 

8  L  l 


(2.24) 


from  which  the  modulus  corresponding  to  the  grating  frequency  used  is  found 
to  be  itk  =  7.  Thus  the  system  of  Fig.  l(as  designed  above, will  convert  the 
21  pulse-position  coded  decimal  tuimbers(Fig.  2.5a)  into  residue  numbers  modulo 


The  next  step  is  to  properly  position  the  output  plane  aperture. 

To  do  this,  we  must  calibrate  the  system.  Since  the  residues  modulo  7  of  the 
decimal  numbers  0,  7  and  14  are  all  0,  we  chose  input  channel  7  for  calibration 
purposes.  With  only  the  input  impulse  on  channel  7  present,  Lho  I'.,  pattern 
in  Fig.  2. lconsists  of  a  dc,  +  1  order,  etc.  spots  of  light.  We  pus i Lion  a 
rectangular  aperture  of  width  7Ax  in  I’2  such  that  the  tic  spot  appears  just  inside 
the  left  edge  of  the  aperture  (when  viewed  from  the  detecLor  plane).  The  -1 
order  spot  is  cut  off  and  lies  just  outside  the  edge  of  this  aperture.  The 
system  is  now  aligned. 

With  the  full  input  plane  1’^  pattern(Fig.  2. 5a) illuminated,  the  output 
plane  1^  pattern  of  Fig.  2.5b  results.  It  contains  2L  vertical  channels  (cor¬ 
responding  to  the  decimal  numbers  0  to  20  on  the  21  input  channels).  On  each 

vertical  channel,  the  horizontal  position  of  the  output  spot  (within  the 

rectangular  aperture  at  F  )  lies  in  one  of  seven  possible  locations.  These 

locations  correspond  to  the  seven  possible  residue  numbers  0  to  6. 

The  residues  modulo  7  of  the  decimal  numbers  0  to  b  are  just  0  to  6 

respectively.  In  these  cases  the  dc  spoL  simply  shifts  across  lho  P  aperture 

from  left  to  right.  For  Liu-  decimal  input  7,  the  dc  spot  now  shifts  out  the 


right-hand  side  of  the  I’.,  aperture  and  the  first-order  spot  enters  Lho 


left-hand  side  of  the  aperture.  Thus  the  correct  residue  number  0  results. 
Use  of  a  high  threshold  SLM  (or  a  thresholded  TV  monitor)  at  can  restore  the 
intensity  of  all  spots  of  light  to  the  same  level. 

To  demonstrate  residue-to-decimal  conversion,  the  system  of  Fig.  2.2.  was 
assembled  with  the  same  lenses  used  in  the  above  experiment.  The  conversion 
performed  was  of  the  residue  number  (1,1,3)  with  moduli  (3,4,5)  into  its 
decimal  equivalent  13.  Only  three  input  channels  are  required  for  this  case. 
The  unit  distance  Ax  for  I’q  and  Pj  of  Fig. 2.2  was  now  0.48  mm.  The  three 
gratings  used  on  the  three  channels  at  were  at  3,4,  and  5  cy/mm  (corres¬ 
ponding  to  the  three  moduli  m^  =  3,  m2  =  4,  and  =  5). 

The  three  input  residue  numbers  were  pulse-position  coded  on  three 

channels  at  of  Fig. 2.2  with  Ax  =  0.48  mm  as  shown  in  Fig.  2\6a.  The  resultant 
output  plane  P2  pattern  in  Fig. 2. 2  is  shown  in  Fig.  2. 6b.  The  horizontal  loca¬ 
tions  of  the  impulses  on  each  output  channel  correspond  to  all  of  the  pos^ifcle 
decimal  numbers  that  could  be  represented  by  the  corresponding  residue  input 
in  the  specific  modulus  chosen.  As  seen  in  Fig.  2.6b,  impulses  from  all  three 
channels  occur  only  at  the  horizontal  location  13.  Thus  the  decimal  equivalent 
of  the  residue  input  (1,1,3)  is  found  to  be  13  and  residue— to— decimal  conversion 
has  been  demonstrated. 
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2.10  DISCUSSION 

The  dynamic  range  of  an  optical  residue  arithmetic  processor  is  the 
product  of  the  moduli  used,  this  is  also  the  largest  number  M  that  can  be 
represented  in  residue.  As  the  number  of  moduli  increases,  each  required 
subcalculation  (recall  that  all  calculations  in  all  moduli  are  performed 

in  parallel)  needs  only  a  dynamic  range  equal  to  the  modulus  chosen.  In  the 
pulse-position  coding  scheme  used,  M  resolution  elements  are  needed  in  the 
input  plane  Pq  of  Fig.  1  and  the  output  plane  of  l'ig.  2  to  represent  the 
full  range  of  decimal  numbers  possible.  Thus  the  system's  most  severe 
bandwidth  and  dynamic  range  requirements  occur  during  the  input  and  output 
conversions,  since  the  full  dynamic  range  of  the  decimal  number  must  be  realized 
in  these  stages. 

The  resolution  Ax  possible  and  needed  in  I’q  of  Fig.  2-1  is  related  to  M  and 
the  input  aperture  A  by 


M  =  A/Ax 


(2.25) 


To  increase  M,  we  must  increase  A  and  decrease  Ax.  The  upper  limit  on  A  is  set 
by  the  lens  system  following  Pj.  The  lower  Limit  for  Ax  is  determined  by  the 
minimum  resolvable  spot  size  and  geometrical  accuracy  with  which  Lho  input 
can  be  recorded.  We  investigated  the  positional  accuracy  of  a  FT  lens  f2.H1 
and  found  it  to  be  better  than  0.06X.  In  general  the  accuracy  and  fidelity  with 
which  input  data  can  be  recorded  limits  the  performance  of  an  optical  processor. 


Thus,  we  expect  Ax  and  positional  inaccuracies  in  recording  the  input  data  to 
limit  the  dynamic  range  of  an  optical  residue  arithmetic  system. 
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The  motivation  for  considering  an  optical  implementation  of  a  residue 
arithmetic  processor  was  to  increase  the  dynamic  range  and  accuracy  obtainable 
in  an  optical  system.  However,  the  approximately  1000  point  linear  SBWP 
presently  obtainable  with  scanners  and  spatial  light  modulators  [2.10]  still  does 
not  permit  more  than  10-bit  dynamic  range  (if  the  pulse-position  coding 
scheme  and  optical  systems  described  earlier  are  used).  Two  approaches  by 
which  increased  input  space-bandwidth  and  hence  system  dynamic  range  can  be 
obtained  are  now  discussed. 

The  first  is  the  use  of  input  numbers  in  binary  rather  than  decimal  nota¬ 
tion.  For  an  input  dynamic  range  of  64,000  we  require  16  bits  to  represent 
the  input  data  in  binary  notation (compared  to  a  linear  input  SBWP  of  64,000 
for  the  pulse-position  coded  decimal  representation  of  the  largest  input 
number).  However,  a  binary/residue  converter  is  now  required.  Such  an  optical 
system  can  be  realized  using  the  joint  transform  correlator  adder  of  Fig.  2.4 
for  B  (equal  to  the  number  of  bits)  cycles.  To  visualize  such  a  system,  recall 
that  N  in  decimal  notation  is  related  to  its  binary  representation  (a^.a^^... 
afi)  by 


B  h 

N  =  T  a.  2b 
n=0  b 

where  each  a.  is  0  or  1.  In  residue  notation  N  modulo  m,  is 
b  i 


|N|  =  (an|2°| 

1  'm,  0‘  'm. 

i  i 


+  a1|21|mi  +  a2|2“| 


.] 


(2.26) 


crom  this,  we  realize  that  N  can  be  represented  modulo  m  as  the  residue  sum 
of  B  numbers  2b  in  residue,  and  that  these  2b  modulo  nu  values  are  known  in 
advance.  In  essence,  we  can  real ize(2. 26) by  adding  0  or  2b  to  a  running 
residue  sum  depending  on  whether  the  corresponding  binary  bit  of  N  is  0  or  1. 
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The  resultant  output  after  the  addition  of  the  most  significant  bit  is  then 
| N | m  in  residue  modulo  nu.  Realization  of  this  conversion  by  the  cyclic 
operation  of  one  adder  of  the  type  shown  in  Fig.  2.4  may  be  preferable  to  a 
cascade  of  B-l  such  adders. 

However,  an  alternate  method  exists  whereby  the  dynamic  range  of  N  can  be 
increased,  without  exceeding  the  assumed  input  linear  SBWP  of  1000. 

This  approach  involves  raster  recording  the  input  data  with  zero-retrace  time 
and  use  of  the  folded  spectrum  optical  spectrum  analyzer  [5.12].  In  the  optical 
residue  arithmetic  version  of  such  a  high  dynamic  range  raster  recorded 
decimal/residue  converter,  pulse  position  coding  of  N  is  still  employed. 

However  K  lines  at  P^  of  Fig.  2.1  are  used  to  represent  N  as 

Nxo  -  “Vo  =  (n-knl)V  <2-27> 


In  other  words,  for  a  maximum  N  =  10,000  and  a  linear  positional  or 

resolutional  =  1000,  we  use  K  =  10  lines  at  Pq  to  represent  N.  If 

0  <  N  <  1000,  a  delta  function  of  light  appears  on  line  K  =  1.  If  1000  <_  N  2000, 

the  delta  function  appears  on  line  K  =  2  at  P_,  etc.  If  we  use  N  =  nm.  +  R 

0  l  nn 

and  assume  KN,  =  rra.  +  R  '  ,  then 
L  1  ra  ^ 


N  -  KN  =  (n-r)m.+R  -R'  , 

L  l  m .  m . 

i  l 


(2.28) 


which  reduces  to 

R  =  (N-KN,)  -  (n-r)m.  (2.29) 

m.  L  1 

1 

when  N,  is  a  multiple  of  m. ,  since  then  R'  =  0. 

L  i  m . 

l 

As  a  numerical  example,  consider  the  conversion  of  N  =  168  to  residue 
modulo  5  assuming  a  maximum  =  50.  In  our  raster  recorded  version  of  the 
input  Pq  pattern  in  Fig. 2.1, the  delta  function  corresponding  to  N  =  168  occurs 
18  spaces  (18Ax)  from  the  left  edge  of  the  k  =  4  line  in  P  .  Since  N  =  50 
is  an  integer  multiple  of  m^  =  5,  N  =  0  to  49  appear  on  line  k  =  1,  N  =  49  to 
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99  appear  on  line  2,  etc.  Then, 


185=  (168-150) j  =  1685-1505  =  1685- 


(2.30) 


The  schematic  diagram  for  such  a  system  is  shown  in  Fig.  2.7.  With  K 

lines  required  at  Pq  to  represent  N,  the  m^  grating  at  is  also  replicated 

for  K  lines  (assuming  1:1  imaging  optics  for  L^) .  A  second  grating  at  a  spatial 

frequency  corresponding  to  the  second  is  recorded  on  the  next  K  lines  at 

(or  alternatively  a  grating  at  one  u  ^  is  present  at  P^  depending  upon  the 

system  scenario).  In  all  cases  only  one  spot  appears  within  the  first  K  lines 

of  the  aperture,  etc.  arvl  thus  as  before  the  horizontal  locations  of  the 

peaks  of  light  at  P_  correspond  to  the  desired  R  .  The  vertical  P_  axis 
2  mi  2 

is  now  partitioned  into  sets  of  K  lines  with  each  set  corresponding  to  one 

R  .  Since  the  two-dimensional  SBWP  of  available  spatial  light  modulators  [  .10] 
i 

is  so  large,  this  utilization  of  available  real  estate  at  P^  is  warranted  if 
the  desired  system  dynamic  range  is  to  be  realized.  One  shortcoming  of  the 
system  of  Fig.  2.7  is  the  low  usable  input  light  level. 

This  can  be  improved  if  the  decimal  input  number  is  used  to  control  an 
oscillator  to  raster  record  the  input  as  a  signal  at  frequency  proportional 
to  N  on  K  lines.  The  2-D  Fourier  transform  of  such  a  raster-recorded  input 
pattern  contains  coarse  and  fine  frequency  axes.  The  location  of  the  output 
peak  along  the  fine  frequency  axis  is  the  desired  residue  modulo  the  input 
horizontal  line  scan  rate [  2. 12 ]  .The  modulus  is  thus  set  by  the  input  line  scan 
rate  and  a  set  of  gratings  is  not  required. 

At  the  output  plane  Pj  of  Figs.  2-1  or  2.7,  another  accuracy  issue  arises. 
Positional  accuracy  represents  no  major  problem  since  only  itk  resolution 
points  per  channel  are  needed.  Ratiier,  the  diffraction  efficiency  of  the  grating 
used  at  P^  must  be  sufficient  to  produce  enough  detectable  light  in  the  M/m^-th 
order  of  the  grating.  The  problem  can  be  relieved  by  use  of  larger  nu  values, 
by  use  of  blazed  or  bleached  gratings,  etc.  The  accuracy  of  the  spatial  frequency 
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of  this  grating  will  also  affect  performance  but  this  is  not  expected  to  be 
a  major  error  source  since  the  grating  need  be  produced  only  once. 

Since  the  aperture  of  the  grating  at  1’^  determines  the  spot  size  in  the 
output  plane  ?2,  we  must  insure  that  this  aperture  remains  above  some  minimum 
value  for  all  beams  emerging  from  different  locations  in  PQ  that  strike  at 
different  angles. 

For  simplicity  in  the  analyses  used,  we  assumed  equal  unit  distances 
AxQ  =  Ax.,  =  Ax  in  Py  and  P,.  In  Fig.  2.1  the  P()  resolution  must  he  superb, 
however  there  is  no  reason  to  retain  fine  resolution  in  P.,  since  only  m.  points 
per  line  are  needed  in  P.; .  Thus,  in  practice  choosing  ne  large  and  ux2  >  AXy 
should  improve  performance.  However,  utilization  of  the  available  space 
bandwidth  (SBU)  and  a  full  system  design  will  determine  the  optimum  values. 

We  demonstrate  by  example  this  latter  point  and  why  the  architecture  of  the 
converters  and  adder  shown  may  be  altered.  Consider  the  input  SBW  require¬ 
ments  for  a  residue  adder  to  add  in  parallel  100  pairs  of  numbers  each  of  which 
is  described  in  residue  by  N  =  11  moduli  (.with  the  largest  modulus  equal  to  31). 
Since  max(m  )  =  31  and  N  =  11,  we  require  a  SI1U  =  343  to  fully  represent  one  number 
in  residue.  To  represent  100  numbers,  we  require  a  SBW  =  34,300.  Thus  to  add 
two  such  pairs  of  100  decimal  numbers,  we  require  an  input  SUM  =  68, 600. 

This  value  may  seem  quite  formidable  and  in  fact  it  is  for  many  systems. 
However,  the  input  SGW  achievable  in  real  time  on  a  spatial  light  modulator [ 2. 10] 
in  a  coherent  optical  processor  easily  exceeds  10^  x  10^  =  10^.  Thus  the  input 
SBW  required  for  the  above  example  leaves  931  ,400  oT  the  1  (/’  possihLc  input 
pixels  unused.  As  this  numerical  example  rather  vividly  demonstrates,  the 
increased  SBW  requirements  of  a  residue  arithmetic  system  represent  no  problem 
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Fig.  2.7  Schematic  diagram  of  a  high  dynamic  range 
decimal  to  residue  converter  using  raster 
recorded  input  data. 
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Fig.  2.8  Large  time  bandwidth  decimal-to-residue  converter. 
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in  an  optical  residue  arithmetic  system.  Thus,  utilization  of  the  large  avail¬ 
able  input  SBW  of  an  optical  processor  will  definitely  require  far  more  elaborate 
input  formats  than  those  described  earlier. 

In  the  numerical  example  presented  above  in  which  N  =  11  moduli  with 

42 

max(m^)  “21,  the  largest  input  number  that  can  be  represented  is  about  2 
Thus,  this  system  with  eleven  moduli  has  the  equivalent  accuracy  and  dynamic 
range  of  a  42-bit  digital  computer. 

The  most  promising  design  we  have  found  for  a  high  dynamic  range  decimal- 
to-residue  converter  is  shown  in  Fig.  2.8.  In  this  system,  the  output  from 
an  LED  source  is  expanded  by  lens  L^  and  passed  through  an  acousto-optic  (AO) 
cell  at  P^.  The  wavefront  leaving  P^  is  then  imaged  onto  a  detector  array  at 
by  lens  Lj. 

To  utilize  this  system  topology  as  a  decimal-to-residue  converter,  we 
temporally  modulate  the  LED  (in  time)  with  a  pulse  of  width  At  occuring  at 
time  NAt,  where  N  is  the  input  decimal  number.  The  light  output  from  the  LED 
and  hence  the  signal  incident  on  the  AO  cell  is  thus 

uQ(t)  =  6 (t-NAt)  .  (2.31) 

We  choose  the  acoustic  time  length  of  the  AO  cell  to  be  m^At  and  describe  the 
signal  in  the  AO  cell  (as  a  function  of  time  t  and  distance  x  across  the  P^ 
aperture)  by 

u^(t)  “  £  6(x-n  M^At+t)  .  (2.32) 

n 

This  represents  a  pulse  of  width  At  and  period  M^At  traveling  in  time  t  and 
distance  x  across  P^. 

The  light  distribution  leaving  P^  and  incident  on  ?2  is  thus  the  product 
of  (2.31)  and  (2,32).  The  detector  array  at  P2  integrates  this  signal  for  a 
time  M^At.  The  output  at  P2  is  thus 
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[  (n+DMj-llAt 

u^(x)  «  £  J  6(t-NAt)  6(x-n  M^At+t)  dt  -  <S[x-(N-n  M^)At]  (2.33) 

n  n  MjAt 

Only  one  pulse  will  be  present  somewhere  in  the  AO  line  when  the  LED  pulses 

on.  Thus,  only  one  of  the  detectors  at  ?2  detects  energy.  Its  location  x 

is  proportional  to  N  is  residue  modulo  as  shown  in  (2.33).  The  integration 

time  of  the  photodetectors  is  M^At  (not  NAt).  The  photodetector  outputs  are 
threshold  detected  and  cleared  each  M^At  seconds.  Thus  detector  noise  which 

will  limit  the  system's  dynamic  range  normally  in  a  time  integrating  correlator, 

does  not  affect  this  system.  The  pulse  itself  present  on  the  detector  is  only 

of  duration  At,  but  the  position  of  the  pulse  within  the  AO  cell  is  not  known 

(hence  the  MjAt  integration  time  at  the  detector).  The  time  aperture  length 

of  the  AO  cell  need  be  only  M^At.  Thus  the  system  is  most  conducive  to  direct 

realization. 

2.11  ALTERNATE  APPROACHES  TO  RESIDUE  ARITHMETIC  PROCESSORS 

The  prior  correlation-based  approach  to  residue  arithmetic  processors  are 
original  and  advantageous  compared  to  others.  However,  for  completeness,  we 
summarize  the  residue  work  of  other  (most  notably  A.  Huang,  et  al.,  [2. 7,2.8]) 
and  an  alternate  mapping  approach  to  residue  computations.  When  two  residue 
numbers  are  operated  on  (added,  subtracted,  multiplied,  or  converted),  a  fixed 
map  or  look  up  table  can  be  constructed  that  automatically  "looks  up"  the  out¬ 
put  value  given  the  inputs  and  the  operation  to  be  performed.  This  look-up- 
table  approach  is  possible  in  residue  arithmetic  because  of  the  limited  dynamic 
range  of  all  operations  (in  any  modulus,  the  accuracy  required  equals  the 
modulus).  Mapping  approaches  are  directly  realizable  using  ROM  (read-only 
memory)  digital  integrated  circuits  as  shown  in  Fig.  2.9.  The  approach  to  a 
larger  system  of  this  type  is  easily  pipelined  to  produce  high  speed  computa¬ 


tions  in  parallel. 
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The  most  difficult  computations  in  any  residue  system  are  the  conversions 
from  decimal  to  residue  at  the  input  and  from  residue  to  decimal  notation  at 
the  output.  In  the  mapping  approach,  the  output  conversion  is  of  paramount 
concern  since  it  cannot  easily  be  pipelined.  This  occurs  because  the  results 
of  one  computation  are  needed  for  the  next  calculation. 

One  optical  residue  arithmetic  system  using  maps  is  shown  [2. 7, 2. 8]  in 
Fig.  2.10.  In  this  system,  the  positions  of  the  horizontal  and  vertical  slits 
at  Pq  and  P^.  select  one  element  of  the  look  up  table  at  ?2  that  is  imaged  onto 
the  output  at  Pj.  This  2-D  map  scheme  is  of  limited  use.  Alterable  1-D  maps 
(Fig.  2.11)  are  preferable.  In  these  systems  in  Fig.  2.11,  we  assume  the 
two  numbers  are  to  be  added.  In  such  schemes  12.7,2.8],  one  of  the  residue  numbers 
(in  one  modulus) to  be  added  is  used  to  alter  the  mapping  topology  as  shown  in  the 
two  examples  in  Fig.  2.11  (adding  1  and  2  to  input  residue  number  modulo  5). 

One  of  the  most  attractive  methods  of  realizing  an  alterable  map  as  in 
Fig.  2.11  is  by  integrated  optical  switches.  An  example  of  such  a  map  to  add 
+3  and  +4  to  an  input  residue  number  modulo  5  is  shown  in  Fig.  2.12.  In  these 
devices,  each  input  is  a  focused  laser  beam  at  an  input  position  proportional  to 
the  Input  residue  number.  The  diagonal  lines  in  Fig.  2.12  represent  reflecting 
surfaces  that  cause  the  location  of  the  output  beam  to  occur  at  a  position 
proportional  to  the  appropriate  output  residue  number  (pulse  position  coded). 

The  system  shown  adds  4  to  the  residue  input  number  modulo  6.  For  example, 
an  input  at  0  exits  at  4  (the  sum  of  0  and  4  is  4) .  An  input  at  4  exits  at 
2  (the  sum  of  4  and  4  is  8  or  2  modulo  6). 

Tsai  at  Carnegie-Mellon  University  [2.13]  has  studied  such  integrated 
optical  switches.  However,  considerable  research  remains  before  such  systems 
appear  realizable.  Issues, such  as  the  cascading  of  such  switches  and  a  large 
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Fig.  2.9  Simplified  concept  of  a  map  using  a  ROM 
(Read  only  memory)  or  similar  element. 
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Fig.  2.10  Optical  system  to  realize  a  fixed  map. 


Fig.  2.11  Simple  alterable  maps  to  add  1  (Fig.  a)  or  2  (Fig.  b) 
to  an  input  residue  number  modulo  5. 
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reduction  in  their  cost  remain  to  be  solved,  before  the  components 
for  such  an  approach  are  available  and  before  a  final  design  can  be  produced. 
Similar  remarks  can  be  advanced  on  other  alterable  map  conceptual  designs 
[2. 7, 2. 8]  using  other  technologies. 

2.12  RESIDUE  MULTIPLICATION 

Addition  and  conversion  are  cyclic  in  a  residue  arithmetic  computer 
(using  maps  as  well  as  our  correlation  approach).  However,  multiplication  is 
not  cyclic  (incrementing  the  input  by  1  does  not  increment  the  output  cyclically 
by  1).  For  example,  in  Fig.  2.13a  we  show  the  map  for  multiplication  of  an 
input  residue  number  by  2  modulo  4.  Residue  multiplication  can  be  made  cyclic 
as  shown  in  Fig.  2.13b  by  performing  an  input  (map  1)  and  output  (map  2)  permu¬ 
tation  leaving  the  multiplication  map  cyclic.  The  general  approach  to  be  used 
for  the  permutations  (maps  1  and  2  in  Fig.  2.13b)  that  yield  a  cyclic  multi¬ 
plication  remains  to  be  fully  analyzed.  Residue  multiplication  can  thus  be 
performed  as  in  Fig.  2.13b  or  by  cascading  M^-l  adders.  Such  a  cascading 
approach  is  feasible  in  a  residue  system  because  the  moduli  are  small  num¬ 
bers.  It  can  also  be  realized  on  a  single  SLM  using  feedback  or  integrated 
optics. 

2.1.3  RESIDUE  ARITHMETIC  MATRIX  MULTIPLIER 

The  basic  residue  arithmetic  element  in  a  phased  array  radar  processor 
is  a  matrix  multiplier.  In  prior  sections,  we  have  described  several  imple¬ 
mentations  of  the  basic  building  blocks  of  a  residue  arithmetic  processor 
(converters,  adders,  and  multipliers)  and  experimental  demonstrations  of 
several  of  these  operations.  From  these  components  the  block  diagram  of  a 
residue  arithmetic  matrix  multiplier  in  Fig.  2.14  follows  directly.  The  pipe¬ 
lined  system  design  shown  in  Fig.  2.14  requires  parallel  computations  in  real 
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Fig.  2.14  Block  diagram  of  a  resiJue  arithmetic  matrix  multiplier 
for  adaptive  phased  arrav  radar. 


Fig.  2. 13  Simplified  diagram  of  a  residue  arithmetic  optical 
matrix  multiplier. 
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time.  Integration  of  the  component  subsystems  using  LSI,  CCD,  Integrated 
Optics,  etc.  is  the  key  to  such  an  approach,  with  the  hope  that  interconnection 
problems  do  not  arise.  Only  added  research  can  answer  such  questions. 

A  second  approach  (Fig.  2.15)  that  we  advocate  utilizes  the  vast  space 
bandwidth  SBW  available  in  an  2-D  optical  processor  and  the  parallel  nature  of 
such  systems.  In  this  case,  the  major  problem  becomes  one  of  how  to  arrange 
and  format  the  input  SLM  data  to  adequately  use  the  available  SBW  product  to 
perform  the  required  matrix  multiplication  with  minimum  complexity  (without 
feedback  and  using  a  single  SLM) . 

2.14  SUMMARY  AND  CONCLUSIONS 

The  need  for  a  novel  and  fresh  processing  approach  to  the  adaptive  phased 
array  problem  merits  the  use  of  a  new  number  system  such  as  residue  arithmetic. 

As  we  have  shown,  residue  arithmetic  offers  the  advantages  of:  parallel  pro¬ 
cessing,  format  control .pipelining,  high  speed,  and  increased  accuracy  and 
dynamic  range.  Thus,  it  has  great  potential  for  use,  in  situations  where  a  high 
computational  load  is  required,  such  as  adaptive  phased  array  radars.  The  major 
problem  areas  in  such  a  processor  are  the  conversions  from  analog  or  digital 
input  data  into  residue  and  back.  We  thus  emphasized  these  problem  areas  in  our 
experimental  work.  We  formulated  a  novel  correlation  approach  that  is  directly 
Implement  able  in  many  available  technologies.  We  have  likewise  reviewed  other 
approaches  using  maps,  integrated  optics,  etc. 

From  this  phase  of  our  program  we  conclude:  (1)  that  residue  arithmetic 
offers  the  fast  parallel  and  accurate  processing  required  for  adaptive  array 
radar;  (2)  that  pulse  position  modulation  should  be  used;  (3)  that  a  correlation- 
based  formulation  is  more  easily  implementable  than  alternate  techniques  (for 
practical  reasons);  (4)  that  added  research  is  needed  on  the  light  losses 
present  in  certain  approaches  and  on  the  input  format  needed  to  utilize  the 
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vast  SBWP  and  parallel  processing  of  an  optical  processor  and  the  associated 
optical  system  to  realize  a  simple  processor;  and  (5)  that  the  state-of-the- 
art  and  expected  time  frame  at  which  cascaded  integrated  optics  switches  will' 
be  available  appears  to  be  far  too  long  for  such  an  approach  to  be  considered 
at  present. 

We  thus  recommend  that  RADC  remain  abreast  of  advances  in  related 
technologies  such  as  SAW,  AO,  CCD,  ROM,  Integrated  Optics,  etc.  and  uses  the 
expected  time  frame  within  which  these  component  technologies  will  have 
matured  to  be  of  use  in  the  fabrication  of  a  residue  arithmetic  and  to  allow 
proper  architectural  designs  for  such  a  system  dedicated  to  the  adaptive 
phased  array  problem.  We  strongly  advocate  that  an  all  optical  approach  be 
pursued  (as  opposed  to  one  in  which  1-D  devices  are  paralleled  and  pipelined) 
in  which  the  input  format  and  SBWP  issues  advanced  in  point  (4)  above  are 
specifically  addressed. 
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CHAPTER  3 

COHERENT  OPTICAL  CORRELATOR  (COC)  FOR  ADAPTIVE  PROCESSING 

3.1  INTRODUCTION 


An  adaptive  array  is  intended  to  cancel  the  effect  of  noise  background 
and  thus  enhance  target  detection  etc.  This  process  involves  two  steps: 
estimation  of  the  ambient  noise  field  and  determining  the  weighting  coeffi¬ 
cients  that  minimize  the  effects  of  the  noise  field.  In  this  chapter,  we 
consider  use  of  a  coherent  optical  correlator  (COC)  processor  to  estimate  the 
noise  field.  We  then  discuss  the  methods  by  which  the  desired  weights  for 
the  array  elements  can  be  obtained  from  this  output  noise  field  distribution. 

In  Sect.  3.2,  we  discuss  the  basic  operation  of  the  COC  system  and  its 
use  in  estimating  the  noise  distribution.  In  Sect.  3.3,  we  introduce  two 
different  COC  architectures  and  describe  the  system  assembled.  Different 
models  for  the  noise  sources,  the  method  used  to  simulate  them  and  the  noise 
estimates  obtained  from  simulations  are  then  presented  in  Sect.  3. A.  In 
Sect.  3.5,  we  consider  the  COC  system  from  a  statistical  view  point  and 
provide  analytical  expressions  for  the  accuaracy  of  the  noise  estimate 
obtained  and  the  available  resolution  of  such  a  system.  In  Sect.  3.6,  several 
important  post  processing  issues  involving  the  estimation  of  the  weights  are 
addressed.  Experimental  results  are  included  in  Sect.  3.7  and  our  conclusions 
are  advanced  in  Sect.  3.8. 

3.2  BASIC  OPERATION 

The  COC  system  which  we  designed  and  fabricated  computes  the  angular 
distribution  of  the  noise  field  as  its  output.  In  this  section,  we  consider 
the  basic  operation  of  such  a  COC  system. 


3.2 


We  consider  an  antenna  system  with  N  elements  that  receive  noise  signals. 

We  assume  the  noise  distribution  to  be  far  from  the  antenna  and  thus  assume 
plane  noise  wave  fronts  incident  on  the  array.  As  a  result,  the  signals  received 
at  different  array  elements  differ  only  by  time  delays.  Since  the  noise  field 
originating  at  different  angular  positions  in  the  far  field  results  in  different 
element-to-element  time  delays  at  the  antenna,  the  desired  N  (0)  angular  noise 
field  distribution  can  be  estimated  by  determining  these  delays  present  in  the 
received  signals. 


Let  us  consider  the  situation  where  the  far  field  noise  consist  of  M 

discrete  noise  sources,  each  at  a  different  angle  0  and  with  a  different  power 

m 

N  .  The  notation  used  is  shown  in  Table  3.1.  We  consider  the  signal  received 
m 

at  the  central  element  of  the  array  as  the  reference  f  (t) .  The  EOC  system 

R 

forms  the  correlation  of  f  with  all  other  received  array  signals  f  .  Using 

K  n 

the  notation  of  Table  3.1,  the  signal  received  at  the  n  th  element  can  be 
related  to  the  reference  signal  by 

£„(t>  ‘  "  C'1** 

m  -  1 

£R(t)  -i;*,  <«•  (3-lb> 

m  «  1 

where 


t  “  k(R-n)d  sin  6  , 
mn  m 


(3.1c) 


and  k  -  2  tt/A.  When  the  reference  signal  is  correlated  with  the  signal  received 
at  the  n  l.h  array  element,  we  obtain 


[?  <oi 

m  m 

MM 


(3.2) 


3.3 


Table  3.1  Phased  Array  Radar  Notation 


N  Number  of  phased  array  elements 

n  A  given  phased  array  element 

f  Received  signal  at  element  n 

f  f  due  to  a  source  at  angle  6 
nm  n  m 

d  Center- to-center  separation  of  array  elements 

X  Wavelength  of  radiation 


k 

M 


m 

e 

m 

s 

m 


Wave  number,  k  «  2w/X 
Number  of  noise  sources 
A  given  noise  source 
Angle  of  noise  source  m 
Noise  source  m  at  angle  6| 


We  assume  that  the  statistical  characteristics  of  the  noise  sources  at 

different  angular  locations  are  independent.  The  cross  correlation  of  f  (t) 

Rm 

and  f  (t)  is  thus  zero  unless  m  -  j  and  (3.2)  simplifies  to 


M 

l  tl 


Rm(t)®  fRm(t_Tmn)] 


M 


l  (t- 


T  )  . 
mn 


(3.3) 


m 

From  (3.3),  we  see  that  the  correlation  of  f  (t)  and  f  (t)  is  the 

k  n 

summation  of  M  noise  signal  auto  correlations,  each  shifted  by  an  amount  r 

mn 

given  by  (3.1).  When  these  N  correlations  are  performed  in  a  multichannel 
optical  correlator,  the  output  plane  light  pattern  of  Fig.  3.1  results.  The 


vertical  axis  in  Fig.  3.1  represents  array  elements  number  n,  whereas  the 
horizontal  axis  represents  the  time  delay. 


Output  plane  pattern  of  multichannel  correlator. 
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If  the  correlation  peak  due  to  each  noise  source  is  sharp  and  if  the  array 

elements  are  uniformly  placed  along  a  line,  the  correlation  peaks  corresponding  to 

noise  source  m  form  a  line  at  an  angle  proportional  to  6_.  To  see  this  we 

c 

evaluate  (3.3)  at  a  fixed  point  xpn  in  the  output  plane  along  the  n  th  channel 
and  obtain 


M 


fR(t)©fn(t)  t.T  (o)  (3.4) 

1  pn  m  r  v 

where  the  contribution  of  terms  other  than  C  (o)  in  the  summation  are  ignored 

P 

because  of  the  sharp  correlations. 


Thus,  Integrating  along  these  correlation  peaks,  we  obtain  an  estimate 
of  the  noise  power  as  a  function  of  angle.  It  should  be  noted  that  a  multi¬ 
channel  correlator  followed  by  an  integration  along  6  provides  information 
only  about  the  noise  power  and  not  the  noise  bandwidth.  Although  N  (6)  does 
not  yield  complete  noise  field  information,  we  often  can  assume  that  the  noise 
bandwidth  is  identical  to  that  of  the  antenna  and  thus  obtain  all  of  the  desired 
information. 


To  provide  the  integration  along  different  6  lines  in  the  output  plane 
of  Fig.  3.1,  three  candidate  methods  are  suggested.  The  first  involves  the  use 
of  a  rotating  slit  whose  light  output  is  Integrated  onto  a  detector  using  a 
2-D  Fourier  transform  lens.  The  disadvantages  of  this  system  are  the  slit 
width  (which  causes  loss  of  resolution)  and  the  mechanical  motion  required 
(which  reduces  the  processor's  throughput  speed).  In  the  demonstration  stages 
of  this  project,  a  slit  is  still  useful. 

The  second  method  uses  a  detector  with  wedged  shaped  elements  with  parallel 
electrical  outputs  mounted  in  the  output  plane  of  the  optical  correlator.  The 


outputs  from  this  detector  array  are  estimates  of  the  light  in  each  sector  and 
hence  of  the  noise  verses  6.  This  system  eliminates  the  mechanical  motion 
disadvantage  of  the  rotating  slit  system,  but  the  number  of  angles  that  can 
be  isolated  is  still  limited  by  detector  technology. 

The  third  method  uses  a  polar  camera.  In  this  system,  the  output  plane 
light  distribution  is  sampled  at  fine  intervals  along  the  two  cartesian 
coordinates  and  the  camera's  supportive  electronic  processor  then  coverts  this 
input  light  distribution  into  (p,  0)  coordinates.  By  integrating  along  p, 
we  can  then  obtain  the  output  noise  estimate  as  a  function  of  0  as  before. 

This  method  may  be  suitable  for  the  final  stage  of  the  COC  system,  but  it 
requires  too  extensive  support  electronics  to  be  included  at  this  point. 

Thus,  for  the  initial  experiments,  a  mechanical  rotating  slit  was  used  and  for 
the  final  system  a  specially  fabricated  wedge-shape  detector  appears  to  be  the 
best  choice. 

A  block  diagram  of  the  COC  system  is  shown  in  Fig.  3.2.  It  is  described 
in  more  detail  in  Sect.  3.3. 

3.3  OPTICAL  SYSTEM 

In  this  section,  we  consider  two  of  the  COC  optical  systems  considered 
and  explain  their  operation.  We  then  describe  the  system  assembled  and  the 
reasons  for  it  as  well  as  the  components  used. 

3.3.1  JOINT  TRANSFORM  CORRELATOR  (JTC) 

The  JTC  is  shown  in  Fig.  3.3.  In  this  system,  the  same  reference  signal 
s  is  recorded  on  all  channels  in  one  half  of  plane  P1  and  the  other  received 
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signals  s^-s^  are  written  in  the  other  half  of  plane  P^  with  a  center-to-center 


separation  of  "2a  .  Since  the  system  uses  1-D  operations,  we  can  consider  what 
occurs  to  each  line  separately.  The  light  transmittance  along  the  n  th  line 


in  P^  is 


Pln(x)  «  sn(x-a)  +  sR(x  +  a). 


(3.5 


Lens  system  L,  forms  the  1-D  Fourier  transform  of  P,  (x)  in  plane  P_.  When 
l  in  / 

this  pattern  Is  recorded  at  P^,  the  subsequent  transmittance  of  Is 

P2n(f>  “  |Sn(f)ej2lrfa  +  SR(f)e'J2lrfa|  ’  (3-6 

where  S  (f)  is  the  Fourier  transform  of  the  signal  s  (x) .  This  recording  is 
placed  in  plane  an^  its  1-D  Fourier  transform  (Ff)  in  x  is  produced  at  Pj 
by  Lj-  At  the  output  plane  P^,  we  obtain 

P3n(f)  “  jR.  {  |  SR(f)  |  2  +  |SR(f)|2}  at  x  -  0 

+  36  (S(f)  S*(f)>  at  x  -  -2a 

n  K 

+  3t{S*(f)  S„(f)}.  at  X  -  +2a  (3.7) 

n  k 


From  (3.7),  we  see  that  the  desired  correlation  is  located  at  x  -  +2a. 
The  three  terms  (3.7)  can  be  separated  in  the  output  plane  if  "a"  is  chosen 
larger  than  2b,  where  b  is  the  width  of  each  signal  recording  in  P^.  The 
disadvantage  of  this  system  is  that  the  full  space  bandwidth  product  (SBW) 
possible  in  P^  is  not  utilized.  Similarly  the  SBW  and  aperture  requirements 
of  the  input  spatial  light  modulator  required  at  P^  and  of  the  lens  system  are 
much  more  stringent  for  this  joint  transform  correlator  topology. 

We  thus  reject  this  system  in  favor  of  the  FPC  system  of  Sect.  3.3.2. 


3.3.2  FREQUENCY  PLANE  CORRELATOR  (FPC) 

The  FPC  system  is  shown  schematically  in  Fig.  3.4.  The  lens  system  and 
1*2  are  1-D  FT  systems.  At  P2»  we  record  on  all  lines  using  holographic 
techniques.  When  the  different  received  signals  S^"SN  are  written  on  different 
channels  in  plane  P^,  the  light  leaving  P2  is 

P2(f)  -  S* (f )  Sn(f).  (3.8) 

Lens  system  performs  a  1-D  FT  of  (3.8)  along  the  f  axis.  This  results 
in  the  desired  multi-channel  correlations  in  the  output  plane  P^.  The  output 
will  contain  other  undesired  terms,  but  these  can  easily  be  separated  from  the 
correlation  by  proper  choice  of  the  signal  to  reference  beam  angle  used  when 
forming  the  P2  pattern.  In  this  system,  the  SBW  of  the  P.^  plane  material  can 
be  much  better  utilized  than  in  the  JTC  system.  This  FPC  system  also  has  less 
stringent  lens  requirements  then  does  the  JTC  one.  Because  of  these  reasons, 
the  FPC  system  was  chosen  for  the  COC  adaptive  processor. 

3.3.3  SYSTEM  DESIGN 

The  lay  out  of  the  FPC  system  that  we  assembled  is  shown  in  Fig.  3.5.  A 
50m  W  helium-neon  laser  (X  «  633nm)  was  used  as  the  coherent  light  source.  This 
light  output  was  divided  into  reference  and  signal  beam  paths  by  the  beam  splitter 
(BS)  and  mirrors  (M)  shown.  The  lenses  used  are  described  in  Table  3.2.  The 
angle  of  the  reference  plane  wave  illumination  was  chosen  to  be  15°. 


Table  3.2 


Component  Specification 


The  P2  pattern  in  Fig.  3.4  is  a  matched  spatial  filter  (MSF) .  It  is 
formed  with  s  recorded  on  all  lines  in  P.  by  recording  the  interference  of  the 
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FT  SR  of  sR  and  the  reference  beam.  For  all  experiments,  the  inputs  at  P^  were 
placed  in  a  liquid  gate  (LG)  and  the  MSF  at  P^  was  recorded  on  a  holographic 
plate.  The  BS  was  orientated  so  that  the  ratio  of  the  reference  power  to 
signal  power  at  the  dominant  carrier  spatial  frequency  was  close  to  unity.  This 
yields  the  maximum  modulation  or  fringe  constrast  in  the  MSF.  A  wedge/ring  detector 
was  placed  at  P2  and  the  BS  was  adjusted  during  MSF  synthesis  to  achieve  equal 
beam  balance  and  modulation  of  the  MSF.  The  exposure  time  was  chosen  to  produce 
a  bias  exposure  with  a  bias  transmittance  of  0.5. 


i 


The  MSF  recording  was  then  developed  in  D-19  solution  for  five  minutes  at 
70°  F.  After  the  emulsion  was  dried,  the  holographic  plate  was  placed  in  a 


hologram  plate  holder  (HPH)  and  the  reference  beam  was  blocked.  The  signals 
s^-Sjj  are  recorded  on  different  lines  in  (on  film  during  initial  experiments) 
which  was  then  placed  in  the  LG.  The  resulting  multichannel  correlation  pattern 
in  n  (3.3)  appears  in  the  output  plane  P^  where  it  is  integrated  in  8  by  a  slit. 
This  integrated  intensity  is  then  recorded  versus  8  by  the  support  electronics. 

3.4  NOISE  MODELS 

In  this  section,  we  describe  the  computer  simulation  used  to  model  the 
COC  system  and  to  model  the  different  noise  sources  and  noise  conditions 
considered  for  the  array  antenna.  The  FORTRAN  routine  that  performs  the 
multi-channel  correlation  is  included  in  the  Appendix. 

3.4.1  SOFTWARE 

The  optical  COC  multi-channel  correlator  was  digitally  simulated  in 
FORTRAN.  The  MAIN  program  generates  single  frequency  jammer  signals  of  the 
desired  amplitude  and  frequency.  The  subroutine  CORREL  performs  the  desired 
correlation  and  PL0T2D  is  used  to  produce  a  2-D  display  of  the  output.  The 
routine  NOISAN  plots  the  noise  distribution  as  a  function  of  angle.  It  utilizes 
the  subroutine  PLOT1D  to  produco  the  1-D  output  displays.  The  subroutine 
COFPEL  produces  the  correlation  of  the  two  sequences  by  forming  the  FFT  of 
each,  the  product  of  these  FTTs,  and  the  inverse  FFT  of  the  result  using  a 
circular  correlation.  By  appending  a  sufficient  number  of  zeros  to  both 
sequences,  the  desired  linear  correlation  results. 

3.4.2  SIMULATION  NOISE  MODEL 

As  our  first  noise  model,  we  consider  the  received  signals  from  one  noise 
source  at  a  single  frequency  and  an  angle  of  60°  to  be  present  at  the  received 


elements  of  an  N  =  39  element  array.  A  total  of  128  samples  of  each  of  the 


Output  of  Fig.  3.6  after  scanning  with  a  slit. 


Time  Delay 
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Fig.  3.8,  is  shown  in  Fig.  3.9.  Its  pattern  is  as  expected  and  displays  two 
prominent  peaks  at  the  correct  locations. 

3.4.4  STOCHASTIC  NOISE  MODEL 

The  signals  we  have  considered  thusfar  are  deterministic  and  do  not  fully 
reflect  the  random  nature  of  the  actual  signals  expected.  In  our  third  simulation, 
we  model  each  of  the  two  noise  sources  by  random  sequences  to  better  simulate  the 
presence  of  two  uncorrelated  noise  sources.  To  produce  such  signals,  a  long 
random  sequence  of  5000  points  was  generated  using  the  RAN  function  in  Fortran. 

The  first  2000  points  were  chosen  as  one  sequence  and  the  last  2000  points  as  the 
other  sequence.  The  resultant  two  signals  are  the  two  uncorrelated  sequences 
desired . 

The  two  digital  random  sequences  obtained  have  white  spectra.  This  was 
confirmed  by  calculation.  The  desired  received  radar  signals  are  band  limited 
noise  patterns.  Since  in  the  optical  experiments,  we  will  record  these  signals 
on  film  using  a  film  recorder  with  constant  modulation  out  to  15  cy/mm  spatial 
frequency,  we  chose  to  model  the  bandpass  noise  signals  with  a  spatial  frequency 
spectrum  between  5-15  cy/mra.  Since  the  film  recorder  spots  are  spaced  at 
intervals  of  10  pm,  the  corresponding  Nyquist  frequency  is  1/20  p  =  50  cy/mm 
and  any  attempt  to  record  higher  frequencies  will  result  in  aliasing.  We  let 
the  digital  radial  frequency  of  +  1.0  it  correspond  to  +  50  cy/mm  spatial  fre¬ 
quency.  To  produce  the  desired  noise  spectrum  with  flat  response  over  the 
5-15  cy/mm  band  and  0  elsewhere,  we  used  a  digital  filter  described  below. 

The  uniformly  distributed  white  noise  sequences  were  passed  through  a  low- 
pass  digital  filter  with  extent  from  -5  cy/mm  (-0.1  n)  to  +5  cv/mro  (+  0.1  tt)  . 

The  low  pass  filtered  output  was  then  multiplied  by  a  sinusoidal  frequency 
10  cy/mm  (0.2  it)  to  shift  the  bandpass  spectrum  to  the  desired  center  spatial 
frequency.  The  design  of  the  digital  filter  is  discussed  below. 
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3.4.5  DESIGN  OF  THE  DIGITAL  LOW-PASS  FILTER 

Following  Oppenheim  and  Schafer,  we  use  a  technique  of  bilinear  transformation 
and  a  Butterworth  filter  to  realize  the  digital  low-pass  filter.  Since  the  cutoff 
frequency  desired  is  0.1  t,  we  introduce  the  condition  that  ripple  in  the  passband 
be  less  than  1  dB,  while  in  the  stop  band  we  require  attenuation  of  more  than 
20  dB,  i.e. 

20  log1Q  |Ha(eJ0,llr)  |  5.-1  (3.10a) 


and 


20  log1(J  |Ha(eJ0*15w)  |  «  -20 


(3.10b) 


The  Butterworth  filter  is  described  by 


|Ha(jn)[2 


l 

l  +  (n/o  )2h 

C 


(3.11) 


where  N  is  the  order  of  the  filter.  Solving  (3.10)  using  (3.11),  we  obtain 
N  =  8  and  fl *  0.3446824.  For  this  value  of  all  specifications  are  met. 
The  poles  of  the  digital  filter  are  found  by  choosing  8  left  half  plane  poles 
in  the  s-plane  that  are  uniformly  spaced  on  a  circle  of  radius  fic. 

To  obtain  H  (z) ,  we  first  find  its  digital  transfer  function  in  terms  of 
the  z-transform  using  the  bilinear  transformation 


„  -  ad-z"1) 

s  -i 

(i+z  n 


(3.12) 


The  resultant  transfer  function  is  very  complex,  but  it  can  be  separated  into 
the  product  of  4  cascaded  functions,  which  are  easier  to  implement  both  from 
software  and  from  hardware  viewpoints.  We  thus  write  H(z)  as 
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H(z)  -  H^z)  •  H2(z)  •  H3(z)  •  H  (z)  (3.13) 

where 

H1(z)  -  0.023615 (1+z-1) 2/ [1-1 .  769z+0.8774z-2] 

H2(z)  -  0.023615(l+z“1)2/[l-1.589z"1+0.6864z-2] 

H3(z)  -  0.023615 (1+z""1)2  /  [l-1.4744z-1+0.5647z-2] 

H^Cz)  -  0 .02 3615 (1+z-1) 2/ [1-1*4 | 89z-1+0. 5057z-2]  (3.14) 

Digital  filtering  is  usually  performed  in  a  computer,  thus  the  white  noise 
sequence  x(n)  is  passed  through  to  produce  an  output 

y (n)  -  1.769  y(n-l)-0.8779  y(n-2)+0.023615[x(n)+2x(n-l)+x(n-2) ]  (3.15) 

We  then  pass  y(n)  through  H,,,  Hj  and  multiply  the  result  by  cos  (0.2  tn)  to 
obtain  the  desired  bandpass  noise  sequence.  This  was  done.  To  verify  the  nature 
of  the  resultant  noise  sequences  obtained,  we  performed  a  2048  point  FFT  on  the 
noise  sequence.  The  magnitude  spectrum  obtained  is  shown  in  Fig.  3.10.  A  noisy 
spectrum  rather  than  a  smooth  one  was  obtained  because  of  the  finite  precision 
of  the  FFT.  However,  we  notice  that  the  spectrum  is  limited  to  the  desired 

frequency  band  and  has  above  20  dB  rejection  in  the  stop  band  and  a  fairly 
uniform  response  in  the  passband . 

The  noise  sequence  as  described  above  was  recorded  on  film  using  our  film 
recorder  and  used  in  the  optical  experiments  described  in  Section  3.7.  We  defer 
further  discussion  of  these  results  until  Section  3.7. 

3.5  ACCURACY  OF  THE  COC  NOISE  ESTIMATE 

In  this  section,  we  analyze  the  COC  noise  estimation  method  from  a  statistical 
viewpoint.  We  define  the  estimation  SWR  when  one  noise  source  is  present  and 
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obtain  the  resolution  performance  when  two  noise  sources  are  present. 

Let  the  receiving  array  consist  of  N  (odd  number)  antenna  elements  and  let 
the  noise  source  be  at  an  angle  6.  If  s  (t)  is  the  signal  received  at  the 
central  element,  the  signal  received  at  the  i  th  element  is 

s1(t)  «  sCt-T^,  (3.16) 

where  i  -  -  (N-l)/2. . .0. . . (N-l)/2  and 

“  i  (d  cos  0)/c.  (3.17) 

When  the  array  elements  are  spaced  on  a  line  at  uniform  intervals  d,  the 
N  optical  correlations  yield 

T/2 

C  (t)  «  /  s  (t)  s  (t  +  t)dt 

1  -T/2  1 

T/2 

«  /  8  (t)  8  (t  +  T  -  T  )dt,  (3.18) 

-T/2  1 

where  T  is  the  integrating  aperture  or  the  time  interval  over  which  the  signal 
is  obtained  (whichever  is  smaller) .  As  an  estimate  of  the  noise  distribution, 
we  integrate  the  output  along  radial  lines  and  obtain  the  estimate 

A  1,1 

N(B)  -  l  C  (Bi) 
i  --N1 

N1  T/2  * 

-  I  J  s(t)  s(t  +  Bi  -  6  i)dt,  (3.19) 

n  — N1  -T/2 

1  * 

*’  and  B  “  (d  cos  6)/c  is  proportional  to  the  correct  angle. 


where  (N  -  1 ) / 2 
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where  we  use  the  assumption  that  Rg(i:)  is  essentially  zero  for  large. values  of 
t.  The  SNR  can  then  be  expressed  as 

T 

SNR  =  NT[  j  (l-|Tj/T)[R^(T)/R^(0)]dT]_1.  (3.24). 

-T 


To  better  understand  the  behavior  of  SNR  with  the  various  parameters, 
we  assume  that  the  signal  has  the  exponential  auto  correlaton  function 


Rg(x)  *=  exp  (-a  |  t  |  )  ,  (3.25) 

where  "a"  is  a  measure  of  the  signal  bandwidth.  Substituting  (3.25)  into  (3.24), 
we  obtain 


NT  r2aT-l+exp(-2aT),-l 
"  T  [ - 4a?T  - 3 


(3.26) 


For  large  time  bandwidth  products,  aT  »  1,  we  can  approximate  (3.26)  by 


SNR  =  NaT  «  (No.  Ant.  Elem.)  (Noise  BW)(Observ.  Time).  (3.27) 

From  (3.27),  we  find  that  the  noise  angle  estimate  can  be  improved  by  either 
increasing  the  number  of  antenna  elements  N  or  by  increasing  the  space  bandwidth 
product  (aT) . 


For  purposes  of  resolution  analysis,  we  postulate  the  presence  of  two  noise 
sources  at  angles  0^  and  We  also  assume  that  these  noise  signals  belong  to 

uncorrelated  random  processes  with  zero  means  and  with  auto  correlation  functions 
R^(x)  and  R2(t) .  The  average  noise  estimate  is  then  given  by 


ri 

N 

-  T  l  j 

[R1Ii(B-6*)]+R2[i(6-6*)]} 

(3.28) 

i  -  -N-1 

where 

* 

6i 

■  (d  sin 

U 

rH 

CD 

(3.29) 

and 

* 

S2 

“  (d  sin 

02)/c. 

(3.30) 

-iv 


si 


If  the  auto  correlation  function  of  both  noise  sources  are  sharp,  the 


average  noise  estimate  will  appear  as  shown  in  Fig.  3.11.  The  Rayleigh 


criteria  can  then  be  used  to  judge  whether  the  two  noise  sources  are  resolvable 


For  the  sake  of  simplicity,  we  assume  6*  «  B*  -  -8^  and  ^(t)  -  R2(t).  Then, 
the  two  noise  sources  at  angular  locations  6  and  -6  will  be  resolvable  if 


N 


i 
i  + 


l  .{R(-i8*)+R(i8*)kl/2 
-N  i 


N1 

l 


{R(0)+R(2iB  )), 


-N1 


i  i  0 


(3.31) 


where  i  °  0  is  deleted  from  the  summation  because  Cq(t)  does  not  provide  any 
information  on  the  location  of  the  noise  source.  Here  we  assumed  that  the  two 
peaks  in  the  noise  estimate  were  resolvable  if  N(0)  was  less  than  1/2  (3dB)  of 
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If  we  assume  that  R(x)  is  sharply  peaked,  we  can  approximate  (3.31)  by 

R(6*)/R(0)  4  N/8 .  (3.32) 

From  (3.32),  the  following  observations  can  be  made: 

(i)  since  R(x)  usually  decreases  as  x  increases,  increasing  the  value  of  N  allows 

★ 

smaller  B  values  to  satisfy  (3.32).  Thus  as  N  increases,  the  sources  can  be 
more  closely  spaced  and  yet  still  be  resolvable. 

(ii)  The  aperture  of  integration  T  does  not  effect  (3.32)  and  as  a  result 

a  larger  time  of  integration  will  not  improve  the  available  system  resolution. 

(iii)  Larger  signal  bandwidths  imply  faster  decay  of  R(B),  thus  allowing 

* 

smaller  values  of  6  to  satisfy  the  condition  in  (3.32). 

3 . 6  POST  PROCESSING 

Thus  far,  we  have  been  concerned  with  estimating  the  angular  distribution 
of  the  far  field  noise.  We  now  consider  an  equally  relevant  aspect  of  adaptive 
array  processing:  the  determination  of  the  accompanying  weights  to  be  applied 
to  the  elements  of  the  phased  array.  In  this  section,  we  consider  two  methods 
to  accomplish  this  and  the  difficulties  associated  with  each. 


We  assume  a  scenario  in  which  there  are  M  noise  sources  s  (t)  and  a  target 

m 

return  s(t)  and  we  assume  that  the  desired  target  is  at  boresight.  The  signal 
received  at  the  n  th  antenna  element  is  thus 


M 

f  (t)  =  s(t)+  l  s(t-x) . 
n  ,  m  mn 

m  “  1 


(3.33) 


When  the  adaptive  array  has  weights  w^,  the  array  output  is 


f( 


t)  =  v  w  f  (t)  -  s ( t )  y  w  +yy  w  s  (t-x  ) 

‘-..inn  L  nLL  nm  mn 

n  *  -N  n  nm 


0  (t)+0  (t), 
s  n 


(3.34) 
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Using  vector  notation,  the  above  condition  can  be  written  as 


[WTCW]U- [WTU] [CW]  =  0 , 


(3.39) 


where  W  is  a  column  vector  [W  „l,  ...  Wx,l3  ,  C  is  the  NxN  matrix  (symmetric) 

T 

with  elements  and  JJ  is  the  unit  column  vector  [1  ...  1]  of  size  N.  The 
optimal  weight  vector  is  then 


^CW  / 

w  =  [-= — 3C  u  -  k'c  u, 
wo 


(3.40) 


where  the  term  inside  the  brackets  does  not  affect  the  SNR  as  can  be  seen  from 
(3.35).  From  (3.40),  we  see  that  determination  of  the  optimum  weights  W  requires 
knowledge  of  the  covariance  matrix  C  and  the  operation  of  inverting  £. 


Case  (ii):  Minimizing  Noise: 

Another  view  point  which  reduces  the  computational  requirements  is  to 

minimize  the  variance  of  0  (t)  given  in  (3.34).  This  is  different  from 

n 

maximizing  the  SNR  as  discussed  above.  In  this  case,  we  assume  that  the 
contribution  0^ ( t)  due  to  the  signal  remains  constant,  i.e. 

Z  w  =  1.  (3.41) 

n 

The  noise  variance  is  then 


Var[f(t)3  =  JJ  Cnjwnwj? 

where 

C  =  Z  R  (t  -t  ) .  ,, 

nj  m  m  mn  mj  (3.36) 

To  minimize  (3.36),  we  differentiate  (3.36)  with  respect  to  each  weight  and  set 
the  derivatives  equal  to  zero.  This  results  in  the  set  of  equations: 


< )-« 
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E  w  =  1  (3.42a) 

n 

Z  w  C  -  0  (3.42b) 

n  n  nj 

for  j  -  -N1,  •••  0  •••  N1. 

Considering  the  (N-l)  equations  in  (3.42b)  and  (3.42a),  we  find  the  optimal 
weight  vector  to  be 

W  =  c”1^,  (3.43) 

T 

where  is  a  modified  matrix  and  =  [1,00  •••0]  .  From  (3.43)  we  notice  that 

we  need  to  compute  only  the  first  column  of  the  inverse  matrix  rather  than  the 

complete  matrix  as  in  case  (i)  in  (3.40).  Since  the  inversion  of  this  matrix 

involves  far  more  computational  load  than  estimating  it,  we  chose  to  reduce  the 

computational  time  by  following  such  a  method.  As  suggested  earlier,  this  method 

will  not  yield  the  best  possible  output  SNR,  but  it  is  more  easily  implemented 
optically. 

We  now  proceed  with  an  approximate  solution  to  (3.40),  The  optimal  weights 
given  by  (3.40)  are  first  rewritten  as 

E  C  ,w  =  k,  for  all  j.  (3.44) 

n  nj  n 

In  the  case  of  a  linear  array  with  uniform  spacing  d  and  single  frequency  noise 
sources,  we  find 


t  -t  *  (2tt/A)  (j-n)d  cos  0 
mn  mj  m 

R  (t)  =  exp  (  j  2x  tv) 
m 

where  v  *=  c/A  and  c  is  the  velocity  of  the  radiation.  In  this  case, 
E  w  E  exp [  j 2irv(2r/A)  (k-n)d  cos  8  )  =  constant  for  all  k. 

n  u  m  m 


(3.45a) 

(3.45b) 

(3.44)  becomes 


(3.46) 
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For  the  case  of  a  single  noise  source,  (3.46)  becomes 


Z  w  exp[j4rr  (vd/X)(k-n)  cos  6  ]  ”  constant, 
n  n  a 


(3.47) 


Ignoring  the  phase  factors  in  (3.47),  we  obtain 


E  exp(-j2irbn)  “  exp( j2irbk) , 


(3.48) 


where 


2rvd(cos  0)/X. 


(3.49) 


From  (3.48)  we  see  that  the  weights  must  be  chosen  so  that  the  Fourier 


transform  (DFT)  of  w  results  in  exp(j2irbk)  for  different  values  of  k.  Once 
n 

a  A 

the  noise  distribution  is  estimated,  8  is  known  and  can  be  used  to  estimate  b 


from  (3.49).  Then,  computing  the  weight  vector  involves  only  an  inverse  DFT 
of  the  right  hand  side  of  (3.48).  This  can  easily  be  performed  in  a  digital 
FFT  post  processor. 


We  should  reiterate  the  above  conclusion  that  a  simple  FFT  post  processor 
is  of  use  in  the  COC  system  only  for  the  case  of  single  frequency  signals  and 
for  uniformly  spaced  linear  array  antennas. 


3.7  EXPERIMENTAL  RESULTS 


In  this  section,  we  consider  the  optical  experimental  computations  of  the 
noise  estimation.  The  data  used  at  and  in  the  system  of  Fig.  3.5  were 
recorded  on  film  using  our  film  recorder.  The  recorded  data  is  described  and 
the  expected  estimates  are  compared  to  the  outputs  obtained  experimentally. 


We  simulated  the  presence  of  two  band  pass  noise  sources  as  described  in 
Sects.  3.4.4  and  3.4.5  and  recorded  these  signals  on  film  in  the  format  described 


<*'*  ' 


j 
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in  Sect.  3.3.  The  two  noise  signals  are  uncorrelated  and  have  a  spectrum  that 
is  limited  to  the  spatial  frequencies  (5-15  cy/mm) .  These  signals  were  recorded 
in  the  format  shown  in  Fig.  3.12. 

The  shaded  part  of  Fig.  3.12  represents  the  opaque  region  of  the  recorded 
film.  The  line-to-line  spacing  used  was  10  pm  and  the  recording  spot  size  was 
10  pm.  The  composite  reference  signal  f^(n,x)  was  recorded  on  1000  lines  in 
the  top  half  of  the  film.  If  the  two  uncorrelated  noise  signals  are  s^(x)  and 
82 (x),  then  f^  is  described  by 

f^(n,x)  “  s^(x)+S2(x)  for  all  n.  (3.50) 

f 

The  received  signals  at  the  other  antenna  elements  are  similarly  described  by 

f2(n,x)  -  s^(x-t-n*10  pm)+S2(x-n-20  pm).  (3.51) 

The  horizontal  locations  of  s^  and  S2  where  chosen  to  correspond  to  a  delay  of 
10  pm  in  one  direction  for  s^  and  a  delay  of  20  pm  in  the  other  direction  for 
82*  The  top  half  of  the  film  is  used  during  preparation  of  the  MSF  at  P2  and 
the  bottom  half  during  the  correlation  process  itself. 

When  f2  is  placed  at  of  Fig.  3.5  and  an  MSF  of  f^  at  P2,  the  expected 
multichannel  optical  output  correlation  plane  pattern  should  consist  of  narrow 
correlation  peaks  at  angles  corresponding  to  the  +  45°  angles  of  the  noise 
sources  in  the  model  as  shown  in  fig.  3.13.  The  central  channel  in  both  films 
was  left  blank  since  it  ioes  not  contain  any  useful  information  and  since  the 
dark  line  present  there  serves  us  well  for  alignment  purposes. 

The  observed  correlation  peaks  from  the  experiment  appear  to  form  into  two 
lines  as  in  Fig.  3.13  at  angles  9^  and  @2  corresponding  to  the  correct  angles 


Fig.  3.12  Format  of  recording  on  the  film. 


Channel 

Number 


-  - V  T 

Time  Delay 


Fig.  3.13  Expected  correlation  output  pattern 
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of  the  noise  sources.  An  isometric  view  of  the  output  obtained  is  shown  in 
Fig.  3.14.  It  clearly  exhibits  the  two  angular  lines  of  correlation  peaks 
predicted  by  theory.  The  results  of  an  angular  slit  integration  of  the  pattern 
in  Fig.  3.14  along  radial  lines  provides  the  N(6)  noise  estimate  shown  in 

Fig.  3.15.  As  shown,  the  output  is  as  expected.  Various  other  cases  yielded 
similar  excellent  agreement  between  theory  and  experiment. 

3.8  CONCLUSIONS 

The  theoretical  analysis  of  a  frequency  plane  correlator  COC  system  for 
adaptive  processing  has  been  performed,  an  optical  system  has  been  fabricated, 
a  digital  simulation  of  the  system  has  been  completed  and  various  noise  source 
models  have  been  produced  and  analyzed.  The  results  of  our  initial  experiments 
using  the  optical  FPC  version  of  the  COC  system  to  compute  the  angular  noise 
field  distribution  and  noise  estimates  were  successful  and  agreed  well  with 
theory  as  shown. 

A  detailed  analysts  of  the  statistics  of  the  noise  estimate  obtained  from 
the  COC  system,  the  postprocessing  requirements,  and  the  number  of  samples 
necessary  to  obtain  an  adequate  estimate  of  the  noise  statistics  were  conducted 
and  provided  the  valuable  information  on  further  direction  for  an  alternate 
COC  processor  for  adaptive  radar. 

From  a  statistical  analysis  of  the  method  of  estimating  the  angular  noise 
distributions,  we  found  that  the  SNR  of  the  estimate  was  directly  proportional 
to  the  number  of  antenna  elements  and  the  space  bandwidth  product  of  the  noise 
signal.  We  also  observe  that  the  resolution  of  the  estimate  was  independent 
of  the  integrating  aperture  and  improved  with  an  increase  in  the  noise  band¬ 
width  or  an  increase  in  the  number  of  antenna  elements. 

We  have  also  determined  that  the  postprocessing  to  calculate  the  adaptive 
weights  from  the  output  of  the  COC  system  can  be  approached  either  by  maximizing 
SNR  (this  requires  knowledge  of  the  complete  covariance  matrix  and  its  complete 


inverse)  or  by  minimizing  the  noise  variance  (this  requires  knowledge  of  the 
full  covariance  matrix  but  only  the  first  column  of  its  inverse).  This  latter 
solution  results  in  a  considerably  simpler  COC  electronic  postprocessor. 

Moreso,  our  analysis  demonstrated  that  the  postprocessor  that  computes  the 
weights  from  the  angular  noise  distribution  outputs  of  the  COC  system  is  a 
simple  inverse  FFT  system  only  for  the  case  of  a  linear  array  and  for  mono¬ 
frequency  noise  sources.  These  conclusions  are  less  disturbing  when  time  and 
frequency  diversity  are  considered,  since  in  this  case,  we  can  simply  bandpass 
the  data  into  narrow  enough  frequency  bands  to  allow  use  of  a  simple  inverse 
Fourier  transform  digital  postprocessing  operation  for  each  adaptive  frequency. 

A  more  major  concern  with  the  frequency  plane  correlator  version  of  the 
COC  adaptive  processor  is  its  limited  input  time  bandwidth  product  and  the 
failure  for  2-D  spatial  light  modulators  to  mature  signigicantly  in  the  past 
year.  Since  2  N  samples  at  the  noise  bandwidth  are  necessary  per  antenn  i 
element  to  accurately  estimate  the  noise  statistics  for  an  N  element  adaptive 
array  and  since  we  must  record  the  heterodyne  carrier  to  obtain  the  necessary 
angular  resolution,  a  very  large  input  time  bandwidth  product  is  necessary 
for  an  optical  COC  adaptive  radar  processor.  We  thus  propose  to  utilize  a 
COC  system  with  1-D  acoustic-optic  transducers  of  high  center  frequency  and 
with  more  reliability  and  availability  than  the  2-D  spatial  light  modulators 
initially  contemplated  and  required  in  the  FPC  version  of  the  COC  system. 

This  will  solve  our  component  availability  problem.  To  obtain  the  large 
time  bandwidth  product  necessary,  we  will  employ  a  time  integrating  correlator 
structure.  We  plan  to  pursue  these  new  COC  system  design  considerations  in 
the  next  phase  of  this  program.  The  basic  conceptual  idea  of  this  new  COC 
system  is  the  same  as  the  one  described  here.  The  modifications  noted  are 
made  to  aLlow  adequate  noise  statistics  to  be  obtained  and  sufficient  angular 
resolution  to  be  retained.  A  version  of  this  time  integrating  acousto-optic 
system  with  both  angle  (space)  and  frequer  y  diversity  also  appears  possible. 
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CHAPTER  4 

ITERATIVE  OPTICAL  PROCESSOR  FOR  ADAPTIVE  PHASED  ARRAY  RADAR 

4.1  INTRODUCTION 

The  most  immediately  realizable  electro-optical  processor  for  APAR  that 
we  have  investigated  is  described  in  this  chapter.  It  is  basically  a  non¬ 
coherent,  feedback,  iterative,  vector-matrix  multiplier  using  a  linear  LED 
input  source  array,  a  fixed  covariance  matrix  and  an  orthogonal  linear  output 
photodetector  array  with  appropriate  intermediate  optics.  The  APAR  loop 
equation  is  reformulated  in  Section  4.2  as  a  vector-matrix  equation.  The  basic 
iterative  optical  processor  (IOP)  system  is  then  described  in  Section  4.3.  The 
system  we  fabricated  to  test  the  components  and  to  demonstrate  use  of  this  novel 
system  design  are  described  in  Section  4.4.  Since  the  data  to  be  processed  are 
complex  numbers,  the  various  methods  to  accommodate  such  operations  in  this 
system  are  summarized  in  Section  4.5. 

Our  experimental  demonstration  of  the  use  of  the  IOP  system  in  an  APAR 
application  on  complex  data  are  then  described  in  Section  4.6.  A  new  electro- 
optical  color-multiplexed  IOP  technique  for  accommodating  complex  data  in  a 
vector-matrix  processor  is  presented  in  Section  4.7  together  with  our  experi¬ 
mental  demonstration  of  this  system.  Several  general  system  application  issues 
are  then  addressed  in  Section  4.8. 

4.2  MATHEMATICAL  FORMULATION 

The  N  intermediate  differential  equations  that  describe  the  average  weights 
W  are  given  by  _ 

IW/G+(M+I/G)W  =  S*.  (4.1) 


■  m 


— - 


The  steady  state  solution  of  Is.  11  for  W  is  obtained  bv  setting  W  "  0  and 


assuming  t. 


1  to  be 


.4  — 1 

MW  =  S*. 


.1  \ 

The  desirod  weights  W  can  bo  obtained  from  S*  and  M  bv 

_  l-» 

W  =  M  S*. 


The  eleetro-oi't  ieal  system  described  in  this  chapter  solves  (u.d)  for  the  weight 
W  using  the  iterative  algorithm 


»  A*  4 

W  =  S  +  tl-Ml  W. 
l+l  t 


where  W.  is  the  ith  estimate  ol  W  (.after  iteration  i  1  .  K.  ,  is  the  i+lst 

l  t+! 

— Sk  — k  A  •k 

estimate,  etc.  When  W.  =  W.  ,,  i-i.-tl  reduces  to  0  -  S  -  MW.  or  S  =  MW. 

t  l+t  t  t 

the  correspond  ing  W  value  is  the  solution  of  t-i.d'l  as  desired. 


s .  .1  1_0_P  SYSTEM  DESCRIPTION 

Refer  to  the  upper  portion  of  the  e  lee  t  ro-opt  ica  1  system  in  Fig.  •< .  1  .  The 

inputs  to  the  system  at  plane  P^  are  the  elements  W^  of  the  vector  W  .  The 

transmit tance  of  the  mask  at  F  ,  is  assumed  for  now  t o  be  M  ,  the  element  s  of 

mn 

the  covariance  matrix  M.  The  svstem  in  Fig.  -..1  is  drawn  to  demonstrate  the 

system's  operat  ion  rather  than  its  detailed  fabrication.  Specifically,  tin- 

input  light  source  at  F.  is  an  arrav  of  l.FFs.  Their  outputs  are  W  .  The  input 

l.F.i''  arrav  is  imaged  vertically  onto  F,  and  the  output  ol  each  IFF  is  expanded 

her  i  cont a  1 1 v  to  uniformly  illuminate  the  corresponding  horizontal  line  ot  the  F 

mask.  The  light  disti  ihut  ion  leaving  F,  is  described  in  .'-F  bv  W  M  .  The 

n  mn 

ont  ical  svstem  between  F  ,  and  F,  (hen  images  !' .  hor  i  cont  a  1  1  v  onto  a  linear  hori- 
‘  .  % 

.rental  phot  odet  oc :  or  arrav  at  F,.  Fach  column  of  F,  is  locused  onto  a  single 
phot  odet  oc  t  or  at  .  The  resultant  light  distribution  at  F,  is  thus 
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Fig.  4.1  Schematic  diagram  of  an  iterative  optical 

processor  for  adaptive  phased  array  radar  processing 
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2  Schematic  diagram  of  the  IOP  svstem  t ahr  iea ted  . 
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We  have  modified  this  basic  electro-optical  vector-matrix  system  first 


described  in  [4.1]  for  use  in  adaptive  phased  array  radar  processing  and  in 

the  solution  of  vector-matrix  equations  by:  (a)  modifying  the  transmittance 

k 

of  P,  to  be  [ I— M] ,  where  I  is  the  identity  matrix,  (b)  addition  of  S  to  the 

output  at  P^  and  (c)  returning  this  output  back  to  the  LED  input  at  P^  [4. 2, 4. 3], 

★ 

The  output  at  P.,  is  now  W  [  I-M  ].  After  addition  of  S  ,  the  output  is 
r  3  n  mn 

k  - 1 

V,'  f  L— M  1  +  S  .  The  corresponding  W  input  that  produced  this  output  :  =  denoted 
n  mil  m 

_ \  _ N 

bv  W..  The  output  from  the  parallel  adder  is  fed  back  to  become  the  new  W 
input  estimate,  which  we  denote  by  .  The  system  shown  in  Fig.  4.1  is  thus 

described  bv 


W.  ,  =  W.  [I-M]  +  S 
i+1  i 


(4.6) 


where  W  denotes  the  LED  system  input  and  hence  the  estimate  of  W  at  Iteration 
i,  whereas  W  denotes  the  output  of  the  system  after  iteration  i  and  the  new 
estimate  of  K  for  iteration  i+1.  As  seen,  (4.6)  agrees  with  (4.4)  and  hence 
the  output  of  the  system  when  W.  =  W  +  c  is  the  desired  adaptive  weights 
solution  in  fa. 3)  to  the  steady  state  equation  (4.2)  version  of  the  adaptive  loop 
described  by  (.4.0. 

4.4  l OP  SYSTEM  DESCRIPTION 

The  1  OP  system  fabricated  used  eight  Howl et t -Packard  LED  inputs  (model 
YDS  2-  •  1 0 1 )  in  a  20.  2  x  2 .  S4  mm"  assemble  lubricated  in-house.  The  M  mask  used 
at  P  ,  lor  this  svstem  had  8  x  8  elements  in  an  4  x  4  mm"  area.  The  width  of 
e.ieii  individual  cel  I  in  the  mask  was  0.3  x  0.3  nsa".  The  10  P  system  fabricated 
is  shewn  sehemat  i  ea  1  iv  in  Fig,.  * .  .  A  photograph  ot  this  svsteni  is  included 
in  r  i  g .  ».  t.  Pul  ;e  width  modul.it  ion  was  used  tor  each  mask  element  ,  with  the 
transparent  mask  area  varied  to  produce  ditterent  desired  gray  levels.  The 
design  o:  the  inte>-r  it  ion  optics  1. ,  -  1.  is  chosen  to  domagnilv  P  vertical  lv 

'  'i  - 
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Fig.  4.3  Photograph  of  the  IOP  system  fabricated. 


by  I/M.J4  ■  1.88  and  magnification  of  P2  horizontally  by  “  1.88.  These  values 

2 

were  set  by  the  7.55  x  A  mm  size  of  eight  of  the  output  detector  elements. 

Referring  to  Fig.  A. 2,  d1  +  d2  »  d3  +  d^.M^  ■  d2/d3  and  M34  “  d4/d3*  The 
focal  length  of  all  cylindrical  lenses  used  was  100  un.  The  imaging  conditions 
1/d^  +  l/d2  =■  l/d3  +  1/d^  -  1/100  require  d^d^/d^d^  -  1.  If  d^  »  d^  and  d2  ■»  d^ 
and  M12  *■  1/M.^,  the  values  d^  -  d^  ■  153  mm  and  d2  -  d^  ■  288  mm  resulted. 

The  amplitude  driver  used  to  amplitude  modulate  all  eight  LEDs  in  parallel 
is  shown  in  Fig.  A. A.  Since  the  output  light  intensity  from  the  Hewlett-Packard 
LEDs  used  is  nearly  linearly  related  to  the  input  current  and  is  relatively 
independent  of  temperature,  the  above  circuit  suffices  for  the  present  experi¬ 
ments.  The  output  detector  circuitry  used  is  shown  in  Fig.  A. 5.  The  uniformity 
with  which  each  row  of  the  P2  mask  is  illuminated  was  improved  to  several  percent 
by  defocusing  L^  slightly.  Use  of  LEDs  with  narrower  cone  angles  can  likewise 
improve  uniformity  of  illumination.  The  fiber  optic  connections  we  plan  between 
P^  and  P2  in  the  next  version  of  this  system  will  greatly  decrease  this  problem. 
Use  of  a  fixed  mask  in  front  of  P2  can  adequately  improve  the  uniformity  of 
illumination  problem.  Cross  talk  between  rows  of  the  mask  at  P2  was  measured 
to  be  as  much  as  A0%.  This  can  be  decreased  by  use  of  guard  rings  on  the  LEDs 
and  mask,  use  of  smaller  LEDs  and  moreso  by  the  planned  fiber  optic  system. 

A. 5  COMPLEX  DATA  HANDLING 

In  practice,  for  the  APAR  problem,  the  vector  and  matrix  masks  used  at  P^ 
and  P2  will  contain  complex  elements.  There  are  several  ways  by  which  complex 
data  can  be  handled  in  such  a  non-coherent  system.  One  of  the  most  attractive 
techniques  appears  to  be  the  representation  of  a  complex  quantity  by  its  three 
projections  [A. A, A. 5]  (at  0°,  120°  and  2A0°)  in  complex  space,  i.e.,  we  represent 
a  complex  vector  element  Ar  by 

A  »  A^  +  A<j>,  +  A£ 

■  A0  exp  0)  +  A^  exp  (j2ir/3)  +  A2  exp  (JAtt/3). 


(A. 7) 


In  this  notation,  the  complex  vector-matrix  product  is  then  described  by 


C  -  B[A] 

-  [A^+'Af  + A*J  [Bg  +  Bj  +  B$] 

-  [AgBg  +"A^B£  + "A£B^]g 
+  [A£B£  +  A~B~  +"A^Bg]^ 

+  [A*B2  +  AJB^  +*A^Bg] j 


(4.8a) 


(4.8b) 


Realization  of  a  vector-matrix  multiplier  and  IOP  system  with  complex  vector  and 
matrix  elements  requires  an  input  array  at  of  length  3M  and  a  mask  at  ?2 
of  size  3  N  x  3  M.  We  describe  such  an  IOP  system  operating  on  complex  data 
for  an  APAR  application  in  Section  4.6.  An  alternate  technique  for  handling 
complex  data  is  described  in  Section  4.7. 


4.6  DEMONSTRATION  OF  THE  IOP  SYSTEM 

To  demonstrate  the  use  of  the  IOP  system  described  in  Section  4.3  on  complex 
data  for  the  APAR  problem,  a  simple  two  element  linear  array  with  one  external 
noise  source  at  an  angle  0^  with  respect  to  the  normal  to  the  array  is  considered 
with  the  addition  of  additive  noise  also  present  [4.2],  The  signals  received 
at  the  two  array  elements  are 

V1  “  X1  +  yl  (4.9a) 

v2  -  x2  +  y2-  (4.9b) 

In  (4.9),  is  the  interference  voltage  in  channel  i  and  y^  is  the  noise  voltage 
in  channel  i.  The  voltage  x2  lags  by  a  phase  angle 

Y  -  (2TTd/A)  sin  8^  (4.10) 

i.e.,  x2  •  x^  exp  (-j y)  and  the  noise  voltages  y^  and  y2  will  be  independent  of 
each  other  and  of  the  interference  signals. 


A.  10 
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We  denote  the  noise  power  In  each  element  channel  due  to  this  external 
noise  source  by  P  and  the  receiver  power  In  each  channel  by  N.  In  this  case, 
the  covariance  matrix  becomes 

[P+N  P  exp  (-JY)”I 

(A. 11) 

P  exp(jy)  P+N  J 

In  the  specific  APAR  problem  considered  to  demonstrate  the  use  of  this  processor, 
we  chose  a  noise  source  emitting  P  -  0.1  watts  at  an  angle  0^  such  that 


(2rd/X)  sin  01  -  tt/3. 


(A. 12) 


We  assume  a  receiver  noise  power  N  «  0.5  (5  times  larger  than  the  received  power 
P  due  to  the  noise  source  itself).  Thus,  N/P  -  5  results  in  a  fast  convergence 
rate  for  the  iterative  processor.  The  actual  mask  used  at  P^  of  Fig.  A.l  is 


[I-M] 


“l-P-N 
P  exp(jAir/3) 


P  exp  (-j  Air/ 3) 
l-P-N 


(A. 13) 


°1’ 


The  form  used  for  this  mask  is  shown  in  Fig.  4.6.  The  matrix  elements  m0', 
m’  are  all  real  and  positive  as  shown  in  Section  A. 5  and  Fig.  A. 6  and  below. 


(4.14a) 

(4.14b) 

(4.14c) 


The  general  form  for  the  mask  is  shown  in  Fig.  A. 6a,  the  specific  form  for  our 
example  is  shown  in  Fig.  A. 6b. 

To  demonstrate  the  use  of  the  electro-optical  IOP  system  of  Fig.  A.l  for  this 
application,  six  of  the  eight  input  LEDs  were  used  and  the  mask  [1  -  M]  was 
arranged  as  shown  in  Fig.  A. 6.  The  six  LED  outputs  at  iteration  I  are  described  by 
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t 


V 


(Wi,0  (1)*  Wi,0  <2>’  Wli  (1>«  wi,l  (2)-  Wi,2  (1)’ 


Wl,2  (2)>* 


(4.15) 


For  our  example,  the  output  from  the  six  photodetector  elements  (after  addition 
of  S)  are  described  by 


(4.16) 


Wi+1  “  (Wi+l,0  Wi+1,0  (2)’  Wi+l,l  (1)’  Wi+l,l  (2)’  Wi+1,2  (1)*  Wi+1,2  (2)) 


This  simple  AFAR  processor  was  simulated.  The  theoretical  results  expected  after 
the  first  six  iterations  in  the  steady  state  (infinite  iterations)  are  shown  in 
Table  4.1.  The  resultant  experimental  results  obtained  (after  the  first,  second 
and  fifth  iterations)  are  shown  in  the  A  scope  displays  in  Fig.  4.7, 

Solving  for  the  adaptive  weights,  we  obtain 


M - 2  2 

L  J  (P+N)  -p 


fP+N-P  exp(-jy)"] 
P+N-P  exp(+jy)J 


We  chose  C  *  0.3,  where  the  steering  vector  is 


(4.17) 


[S 


(4.18) 


This  choice  of  C  simplifies  the  multiplicative  factor  in  (4.17)  to  be  approximately 
unity.  Evaluation  of  the  six  components  of  in  (4.17)  for  the  first  six 

iterations  are  given  in  Table  4.1.  The  evaluation  of  (4.17)  for  this  case  yields 


rui  „  OJL  To. 6-0.1  exp(-jTr/3)-J 
1  J  0.35  (.0. 6-0.1  exp(+jir/3)J 


(4.19) 


The  results  shown  in  Fig.  4.7  are  quite  close  to  the  predicted  values  in  Table  4.1 
as  seen.  As  demonstrated,  the  results  of  this  initial  IOP  system  demonstration 
are  most  promising. 


idetec  t 
i,  1,  a 


4 . 7  COMPLEX  DATA  REPRESENTATION  BY  COLOR  MULTIPLEXING 


The  complex  data  representation  scheme  noted  in  Section  4.5  and  experimentally 
demonstrated  in  Section  4.6  requires  three  times  the  space  bandwidth  product  for 
the  input  LED  array  and  the  output  photodetector  array  and  nine  times  the  space 
bandwidth  for  the  matrix  mask  at  These  input  and  output  plane  requirements 

can  be  reduced  by  use  of  a  novel  color  multiplexing  technique  shown  in  Fig.  4.8. 

In  this  method,  three  linear  input  LED  arrays  (each  of  a  different  color  or  wave¬ 
length)  are  used  and  three  linear  output  photodetector  arrays  are  required.  Each 
of  the  three  components!^  ,  W«  ,  W*  of  H  are  represented  by  a  separate  linear 
LED  input  array.  A  grating  inserted  between  P 2  and  enable  separation  of  the 
three  output  products  on  the  three  linear  phatodetector  arrays.  Summation  of  the 
proper  photodiode  outputs  at  P^  yields  the  correct  vector-matrix  product  which 
can  then  be  fed  back  to  P^  in  the  IOP  system  [4.3]. 

The  results  of  our  initial  experiments  to  demonstrate  [4.3]  the  use  of  this 
color  multiplexed  IOP  system  of  Fig.  4.8  are  shown  in  Fig.  4.9.  In  this  experi¬ 
ment,  an  arc  lamp  source  was  used  with  three  spectral  lines  (A^  violet,  A^  green, 
and  Aj  orange) .  A  two  element  input  vector  with  components 

au  =  (1,0),  aX2  =  (1,1),  aX3  =  (0,1)  (4.20) 

was  used  as  the  2x3  input  (Fig.  4.9a).  We  consider  the  multiplication  of  A  by 
B.  The  output  from  the  2x6  mask  B  used  is  shown  in  Fig.  4.9b.  The  three  rows 
of  the  matrix  B  are  read  out  with  A^  +  A^,  A2>  and  A^  light,  respectively, 
because  of  the  input  pattern  chosen  in  (4.20).  The  mask  P2  output  shown  are  as 
expected.  The  P^  output  consists  of  three  rows  of  six  detector  elements  per  row 
(Fig.  4.9c).  These  correspond  to  the  eighteen  possible  cross  products.  All  are 
of  the  proper  form  as  indicated.  These  Initial  experiments  using  this  novel 
complex  data  representation  method  were  most  promising.  Their  future  use  requires 
the  need  for  several  parallel  arrays  of  laser  diodes,  which  we  hope  to  have  and 


Fig.  4.8  Schematic  diagram  of  a  color-multiplexed  optical  processor. 


Kig.  4.9  Color  multiplexed  complex  vector-matrix  multiplication. 

(a)  input,  (h)  mask  output,  (c)  system  output.  (Printed 
in  black  and  white.  V  =  violet,  (’.  =  green,  0  =  orange). 
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demonstrate  in  Phase  2  of  this  program. 

4.8  SYSTEM  CONSIDERATIONS 

Either  LEDs  or  laser  diodes  can  be  used  as  the  light  source.  Laser  diodes 
are  preferable  because  of  their  smaller  size,  higher  power  output,  variety  of 
wavelengths  available,  and  their  large  linear  dynamic  range  of  output  light 
intensity  with  input  current.  Laser  diodes  require  heat  sinking  to  compensate 
for  temperature  effects,  but  often  this  can  be  performed  directly  on  the  chip. 

Pulse  width  source  modulation  appears  to  be  necessary  with  LED  sources  because 
of  their  non-linear  behavior  at  low  current.  This  will  decrease  the  system's 
speed  but  it  is  still  more  than  adequate  for  the  APAR  application.  A  system  cycle 
time  of  1  psec  appears  possible.  However,  multiplexed  input  and  output  electronics 
are  preferable  for  reasons  of  cost  and  man  hours.  Even  with  pulse  width  modula¬ 
tion  and  input  and  output  multiplexing,  adequate  system  speeds  for  the  APAR 
program  are  easily  obtainable. 

A  major  concern  in  such  an  iterative  processor  is  its  accuracy.  The  results 
of  computer  simulations  have  shown  that  cross  talk  and  nonuniform  illumination  of 
the  P2  mask  and  detector  nonuniformity  are  the  major  error  sources.  As  suggested 
in  Section  4.4,  use  of  fiber  optic  connectors  between  P^  and  P^  appears  to  solve 
the  uniformity  and  cross  talk  problems.  A  fixed  auxiliary  mask  in  front  of  P2 
is  planned  to  compensate  for  any  nonuniformity  that  still  remains.  The  optics 
between  P2  and  will  be  removed  in  the  next  system  by  designing  the  front  end 
so  that  the  detector  array  can  be  mounted  flush  with  the  P2  mask.  LED  nonlinearity 
will  be  decreased  by  use  of  pulse  width  input  modulation.  Any  differences  that 
remain  between  the  saturation  levels  of  the  different  LEDs  and  any  sensitivity 
differences  present  between  individual  photodetectors  will  be  compensated  for 
electronically  using  a  programmable  read  only  memory  in  the  electronic  feed¬ 


back  circuitry. 
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This  IOP  system  appears  to  be  most  promising  and  to  merit  further  research 
with  major  attention  given  to  its  accuracy  and  to  fabrication  of  the  prototype 
system  utilizing  recent  advances  in  LED,  fiber  optics,  detector,  and  solid-state 
electronics  and  the  extensive  test,  analysis  and  demonstration  of  the  system 
fabricated . 
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CHAPTER  5 


SUMMARY  AND  CONCLUSION 


Three  different  and  new  optical  signal  processing  techniques  and  system 
architectures  for  adaptive  radar  processing  have  been  described,  analyzed, 
initial  experiments  performed,  component  availability  assessed  and  recommenda¬ 
tions  provided. 

5.1  OPTICAL  RESIDUE  ARITHMETIC 

The  first  technique  considered  utilized  residue  arithmetic.  The  two  moti¬ 
vations  for  considering  such  an  approach  for  adaptive  radar  processing  were  the 
high-speed  and  parallel  nature  of  a  residue  arithmetic  processor  (specifically, 
no  carries  are  necessary)  and  that  high  computational  accuracy  is  possible  with 
modest  accuracy  in  each  subcomputation.  Application  of  residue  arithmetic 
techniques  to  an  optical  processor  is  thus  attractive  because  of  the  parallel 
nature  of  an  optical  system,  which  provides  a  vehicle  with  which  to  perform  the 
many  parallel  computations  necessary  in  a  residue  arithmetic  system.  It  is  also 
attractive  because  all  individual  computations  need  be  performed  to  an  accuracy 
equal  to  that  of  the  basis  rather  than  the  full  number.  Thus,  an  optical  residue 
arithmetic  processor  provides  an  attractive  method  to  improve  the  accuracy  of  an 
analog  processor. 

We  first  developed  a  new  correlation  formulation  of  residue  arithmetic 
processors.  This  allows  such  processors  to  be  realized  using  correlators  rather 
than  conventional  systems  that  perform  additions,  subtractions,  etc.  Since 
optical,  CCD  and  other  technologies  perform  correlations  rather  than  the  more 
conventional  mathematical  operations,  such  a  formulation  leads  more  directly  to 
implementation  of  residue  processor  architectures  in  such  technologies. 

As  the  first  two  special  components  of  an  optical  residue  arithmetic  processor, 
we  considered  the  input  and  output  problems,  specifically  converting  digital  and 
decimal  data  into  residue  form  for  input  to  a  residue  arithmetic  processor  and  then 
techniques  for  converting  residue  data  into  digital  and  decimal  form  at  the  output 
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of  the  system.  Optical  decimal/residue  and  residue/decimal  converters  were 
designed,  and  experimental  laboratory  system  was  assembled  and  experimental 
demonstration  of  these  two  new  systems  were  obtained.  An  optical  residue 
arithmetic  adder  and  vector-matrix  multiplier  system  were  then  designed.  ' 

The  two  major  issues  of  concern  in  such  a  processor  are  how  to  input  data 
to  such  a  system  in  the  proper  format  and  the  availability  of  the  necessary 
component  technology.  The  first  issue  requires  more  extensive  study  and 
analysis  as  well  as  a  reformulation  of  the  adaptive  radar  problem.  In  all  sys¬ 
tems  considered,  the  control  system  to  provide  input  data  of  the  proper  form 
and  in  the  proper  locations  appears  to  require  such  extensive  logic  that  the 
parallel  processing  features  of  a  residue  arithmetic  processor  become  less 
attractive.  Implementation  with  2-D  spatial  light  modulators,  integrated  optical 
systems  and  other  optical  technologies  appears  to  require  such  extensive  device 
and  component  development  that  implementation  of  such  systems  appears  to  be  years 
away.  The  most  readily  implementable  residue  arithmetic  processing  technique 
appears  to  utilize  fiber  optics  and  thus  capitalize  on  rapid  advances  in  this 
area  for  communication  applications.  Because  the  other  two  optical  processors 
we  considered  for  adaptive  radar  applications  can  be  realized  with  existing 
components  and  appear  to  be  capable  of  the  necessary  system  performance,  we 
recommend  that  one  keep  abreast  of  developments  In  integrated  optics,  CCD  devices, 
Josephson  junction  technology  and  other  candidate  residue  arithmetic  system 
components  and  that  further  research  be  performed  on  the  COC  and  IOP  systems 
for  adaptive  phased  array  radar. 

5 . 2  COC  PROCESSOR 

The  second  candidate  optical  system  considered  is  the  closest  of  any  to 
a  conventional  optical  processor.  The  basic  system  is  a  hybrid  processor  con¬ 
sisting  of  a  frequency  plane  correlator  modified  to  function  as  a  multi-channel 
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system.  The  multi-channel  output  correlation  plane  pattern  is  then  slit  inte¬ 
grated  using  a  specially  shaped  solid-state  detector.  A  simple  digital  Inverse 
DFT  postprocessor  can  then  compute  the  necessary  adaptive  weights  from  this 
output  for  many  scenarios. 

This  COC  processor  operates  on  the  heterodyned  received  signals  at  the 

antenna  elements.  It  computes  the  correlation  of  the  N  received  antenna  signals 

with  a  reference  signal  from  one  antenna  element  in  parallel.  The  multi-channel 

optically  computed  output  correlation  plane  pattern  is  such  that  each  noise  source 

at  a  different  angle  produces  a  set  of  correlation  peaks  that  align  in  one  row  at 

a  given  angle  in  the  optical  correlation  plane  proportional  to  the  azimuth  angle 

0  of  the  noise  source.  Thus,  after  integrating  the  multi-channel  correlation 
ra 

output  along  different  radial  lines  (using  a  special  solid-state  detector  with 

angular  shaped  detector  elements) ,  the  angular  distribution  of  the  noise 

field  N  (0  )  is  obtained, 
m  m 

For  a  single  monofrequency  noise  source  and  a  linear  array  with  equally 
spaced  elements,  an  inverse  DFT  postprocessor  has  been  found  to  be  sufficient  to 
com;  :**•  the  weights  from  the  N  (0)  information.  Initial  experiments  and  computer 
simulations  of  this  system  produced  outputs  as  expected.  A  statistical  analysis 
of  this  system  was  conducted  in  which  the  accuracy  of  the  estimate  was  related 
to  the  number  of  adaptive  elements,  signal  bandwidth  and  observation  time.  Compu¬ 
tation  of  the  optimum  weights  from  the  N  (0)  pattern  were  performed  using  maxi- 
zation  of  the  output  SNR  and  minimization  of  the  output  noise  as  two  goodness 
criteria.  To  obtain  the  optimum  weights  using  SNR  as  the  criteria  required 
computation  of  the  full  covariance  correlation  matrix  and  its  full  inversion. 
Conversely,  the  second  criteria  required  computation  of  only  one  column  of  the 
correlation  matrix,  as  the  COC  system  performs.  Thus,  we  expect  the  weights 
computed  from  the  hybrid  optical/digital  COC  system  to  be  optimum  for  minimizing 
the  output  noise  rather  than  the  output  SNR. 
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Our  major  concern  with  Che  present  COC  system  was  the  need  for  real-time 
and  reusable  2-D  spatial  light  modulators  and  the  need  to  process  long  time 
bandwidth  product  signals  to  obtain  adequate  noise  statistics.  We  thus  advocate 
use  of  the  same  basic  correlation  and  system  concept  using  instead  a  time 
integrating  acousto-optic  correlator  to  achieve  the  long  signal  integration 
time  bandwidth  products  needed  and  to  utilize  more  available  real-time  trans¬ 
ducers  . 

5.3  IOP  SYSTEM 

The  third  optical  system  for  adaptive  phased  array  radar  processing  appears 
to  be  the  most  attractive  one  that  is  both  most  directly  analyzed  and  most  readily 
implemented.  The  basic  system  uses  a  linear  input  LED  array,  2-D  mask  and  linear 
detector  array  to  perform  a  vector-matrix  multiplication.  Addition  of  a  vector 
to  this  system's  output  and  returning  the  output  to  the  LED  inputs  results  in 
an  iterative  optical  processor  using  feedback.  This  IOP  system  can  then  solve 
the  vector-matrix  equation  of  adaptive  radar  for  the  output  weights  using  an 
iterative  algorithm.  The  inputs  to  this  non-coherent  optical  processor  are  the 
covariance  matrix  M  and  the  steering  vector  S. 

Three  techniques  that  enable  this  system  to  operate  on  complex  data  have 
been  developed  and  experimentally  demonstrated.  These  include  a  new  color 
multiplexed  technique. 

The  accuracy  of  this  or  any  such  feedback  or  iterative  system  was  our  primary 
concern.  Uniform  illumination  of  each  row  of  the  mask,  crosstalk  between  rows, 
and  non-linear  LED  responses  were  determined  by  analysis,  experiment,  and 
simulation  to  be  the  major  error  sources  of  concern.  The  use  of  fiber  optic 
coupling  from  the  LEDs  to  the  mask  (and  a  fixed  mask  in  front  of  the  variable 
mask)  overcomes  the  nonuniform  illumination  and  crosstalk  problems  noted  above. 

The  use  of  pulse  width  modulation  of  the  LED  source  outputs  or  amplitude  modulation 
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of  laser  diode  sources  overcomes  Che  source  nonlinearity  pro! lems  noted . 
Differences  in  the  saturation  level  of  the  LED  sources  and  differences  in 
the  sensitivity  of  the  output  detector  elements  can  be  corrected  in  the  elec¬ 
tronic  feedback  system  itself. 

The  system  fabricated,  initial  tests  performed  on  it,  and  the  results  of 
the  simple  adaptive  phased  array  radar  problems  performed  to  demonstrate  the 
use  of  the  system  were  most  encouraging.  We  thus  recommend  the  fabrication 
of  a  prototype  of  such  a  system  as  described  above  with  extensive  component 
tests  and  attention  to  the  accuracy  and  system  performance  obtainable  together 
with  extensive  testing  of  the  use  of  the  system  and  diverse  adaptive  phased 
array  radar  problems. 
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0010U  C 
0020U  C 
00300  C 
00400  C 
00500  C 
00600  C 
00200  C 

ooeoo  c 

00900  C 
01000 
01100 
01200 


02400 
02500 
02600 
02700 
02800 
02900  C 
03000  C 
03100  C 
03200  C 
03300  C 
03400  C 
03500  C 
03600 
03700  C 
03800 
03900 
04000 
04100 
04200 
04300 
04400  C 
04500  C 
04600  100 

04700 
04800 
04900 
05000 
05100 
05200 
05300 
05400 
05500 
05600 

05700  150 

05800 
05900 
06000 


rwiiinAM  gives  tip-  ukhplat!hn  ok  thk  noise  receiver  for 

Kri!  F NT  nK  AN  ANTf  : M A  AsoAV  v I TH  T"E  N^IGE  RECEIVER  POP 
THK  ELEMENT  IN  TIIK  MIOP.K  UF  TIIK  ARRAY  WHICH  l  WILL  CALL 
PKFEKFNCK  ELEMENT/  THF  1UTPUT  IS  IN  A  PLOTTING  FORK. 

TO  EXFr»»TP  IT  T«E  FOLI.C^ING  SEOHENCF  I'AS  TO  BF  UONF: 

0  TAKE  SYS:TE*6 

0  E*E  !JANTENN.CMr,SVS:TFKT.PEL 

IN  WHICH  ANTENN.CMD  IS  A  COMAND  FILE  WITH  ROB- FOP/ SUB. FOR, 

IM ArF. .FOt?  AND  PLOT. FOB  IN  IT. 


ORAL  At  Lf  011  771/ BEF(0:1  27),  APRAYt  0:127/ 0:38  ),tIPPFPv,t)PPEPN 


01300 

01400 

INTEGEP 

A/X/N 

01500 

A  = 

1 

01600 

upopoy  = 

177.0 

01700 

nppERN  = 

38.0 

01600 

FBI  = 

0.1 

01900 

FP2 

0.09 

02000 

AN^LEl  = 

0.5 

02100 

ANGLE2  = 

0.707 

02200 

MAX  = 

0.0 

02300 

PI 

4* ATAN { 1.0) 

IN  THIS  PROGRAM  1  GAVE  TWO  POSSIBLE  INPUTS  AS  NOISE  SOURCES 
AT  FIPST  THE  NFXT  00  LOOPS  IN  THE  PROGRAM  CAN  BE  USED  TO 
£«nt.ATE  T«E  HHIFE  0®  A  SOURCE  W*»Iph  FIVE  A  WAVE  AS  rOS(PHI) 
ONLY  FOB  phi  BETWEEN  ZERO  AND  TWO  PI/  IT  MEANS  ONLY  ONE  PERIOD 
ITS  A  NON-PcRTCDICAL  NOISE. 

THE  OTHrR  POSSIBILITY  IS  A  SECOND  CPO"P  r>0  LOOPS  FOP  THE 
CASE  WHEN  THE  NOISE  IS  PERIODICAL(COS(PHI ) ) . 


NON-PERIODICAL  NOISE 
no  ICO  X=0/ 

REF( X) =0. 0 
H° 1  =  2*FT*FR1*X 
HR2  =  2*PI*FR2*X 
TF( (HR1.GE. 
I  °E  F ( X  )  = 

I F { (HR7.GE. 
1  A  * ( COS( 

CONTINUE 

PO  200  N=0, 

PO  l'O 

ALL( 
H  A 1  = 

HA2  = 


IJPPFRX 

-  FI*ANGLE1*19 

-  PI«ANGLE2*1P 

0.0 I. AND. (HP1. LE. ( 2*PI ) ) ) 
A*(CQS((2*PI*FR1*X)*CPI*ANGLE1*19))*1) 
0.0).AND.(HP2.LF.(2*PI)))REP(X)-=REF(X)  + 
(2*PI*FR2*X)+fPI*ANGLE2*19))+l) 

UPPCRK 

X=0/ UPPERX 
X)=0.0 

2*PI*FP1*X  -  FI* ANGLF1*N 
7*PI*FR2*X  -  PI*ANGLF,2*(UFPF.RN-N) 

TF( (HAl.GE.O.O).AND. (HA1.LE.(2*PI)))ALL(X) 
A*(COSC(2*PI*FPl*X)*(FI*AKCLEl*N))*l) 
IF((HA2.CE.O.O). AND. ( H A  2. LE . ( 2*P I ) ) ) ALL( X ) 
ALI.(X)*A*(COS(  f2*PI*FR2*X)«' 
(PI*ANGLE2*UFPFRN-N))*1 ) 


CONTINUE 


'>200  C 
u>  300  C 

6400  C 
>500  C 
o600  C1000 

'6100  C 
6800  C 
'6900  C 
17000  C 
■7100  C1500 

7200 
7300 
7400 
7500 
7600 

/  700  200 

7800  C 
7  900  C 
HOOO  C 
'100  250 

8200  C 
9300 
9  400 
9500 
.icOO 
9700 

:eoo 

3900 
7000 
>100 
0200 
,'JOO 
9  400 
•  aOO 
9600 
#700 
9800 

9900  300 

uOOO 
0100 
j200 

030( 

0  40C 
0500 
9600 
>700 

0750  410 

udOO  400 
0900 
1000 


PFRinni'-AT.  NDTS5, 

DU  1000  X  =  0,  UPPER  '• 

!»FF(X)  =  AMCnS(2*P'»FPl*X  ♦  P I  *  ANGLE  1  *19)  ♦ 

1  CQS(2*PI »l K?*X  «■  PI*ANGLB7*19)  ♦  2) 

HINTIN'*5: 

DU  200  N=0,UPPERN 

DO  1500  X=0, UPPFRX 

ALL(X)=A*(CnS(2*PI*FRl *X  ♦  n  I  *  ANCf.F  1  *N  )  + 

1  COS  (  2  *PI  *FR2*X  ♦  PI*Af!GLE2*(U0FEPM-N))42) 

CONTI N"E 


CALL  CORRET (127,ALL,°SF) 

DO  200  X=0, UPPERX 

ARCAV(y/N)=ALL(X) 

CONTINUE 

CALL  PLOT IN(127, 39, ARRAY, 'RORERT.DAT', MAX) 

OPFN ( UNIT  =2 1 ) 

VRITFC>l,?rO)APRAV 
FORMAT!  9UXG10.5)) 

CLOSFCi'NIT-21  ) 

CALL  PLOT  2D ( 127/38/ARRAY) 

PA"SF 

CALL  NOISAN(UPPERX,UPPERN, ARRAY) 

STOP 

END 

T  K  i  \  i  *  nuOTMi 

SUBROUTINE  NOISAN(UPPERX,UPFERN,ARPAY)  , 

REAL  °AC/ATIT  A  (180), ARRAY  (0:127,0:39)  d  -\U  CcmiAA  uoiS 

INTEGER  X0,Y0,X,TITA,  Y,Y1  j'  0 

Pl=  4*ATAN(1.0) 

XO-  uPrER  X /2 
Y0=UPPFRN/2+4 

no  300  TITA-1, 180 

ATITA(TITA)=0.0 

CONTINUE 

no  400  Tir5=l,lP0 

RAD=(TITA/180.0)*PI 
Ic  (TITA.EQ.90)  GCTC  400 
00  410  X=0,U°PERX 

Y=(SIN(RAD)/COS(RAD))*(X-XO)+YO+C.5 

Yl=T!FPERN-Y 

IF((Yl.LT.O).QR.(Y.LT.O))  GOTO  410 
ATITA(TITA)=ATITA(TITA)*AORAY(X,Y) 

CONTINUE 

CONTINUE 

CALL  °f.r|Tlr'(1  "O,  AT  it  A,  I  .  0,  1  .0) 

STOP 
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00100 

C 

00200 

c 

00300 

SUilPUUTINE  CORPELCNU  ,TX,TY) 

00400 

c 

F REFORM*!  A  CIRCULAR  'TRRELATICN  CF  THF  RF.AL 

00S00 

c 

ARRAYS  TX(0:  NUM)  TY(":MUM).  THF  RESULT  IS  I II 

00600 

c 

Ty(0:NU'‘)  A*T  IR  RRIFTF.n  SO  T*1  AT  TX(0  )  =  TX(N'iv/2) 

00700 

c 

OOEOO 

INTEGER  I/MUR/N 

00900 

RFAL  TX(0: 1023 ),TY(0: 1023) 

01000 

Cny°(.EX  X(0:1023),Y(0:1023) 

01100 

DU  400  1=0, MUM 

01200 

400 

X(I  )=CMPLX (T  X (I ), 0. 0 ) 

01300 

no  *00  1=0, vow 

01400 

500 

Y(1)=CMPLX(TY(I),0.0) 

01500 

CALL  FFT(UUy,X,l) 

01600 

CALL  FFT (NUM,Y, 1 ) 

01700 

CALL  ROVJfNUw,Y) 

01800 

DU  600  1=0, NUM 

01900 

6uC 

X( I  )  =  X(  I  )*Y(  I  )/(NUM**2 ) 

02000 

CALL  FFT (NUM,X,-1 ) 

02100 

N  =  N'*M/2 

02200 

DO  700  1=0, N 

02300 

TX(I)=APS(T(T*N*1)) 

02400 

700 

TX(I+N+1)=ABS(X(I>) 

02500 

flFT»'RN 

02o00 

END 

-1 
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00100 

C 

00200 

C 

00300 

SUHPQUTI  KF  FFT(NUM/>  1  N V > 

00  400 

C 

tuis  ^HnpnuT i mf.  rrFr  t,,f  fft  pf  tpe  one 

niv^NnlPNAL 

00500 

C 

CUMf’LE  X  VECTOR  v  (  0  :  N 'JM  ) ,  NUM  =  2  *  *N - 1  .  THIS 

I  S 

00600 

C 

A  rFC!MATT0N-IN-TIf'c/T5Ar'tX  2,IN-PI.ArE  FPT 

ROUTINE 

00700 

C 

TU  HU  INV°SSE  FFT  SFT  1NV  TO  -1  OTHERWISE 

SFT 

TO  1 

00800 

C 

00900 

INTEGER  NUM/ I / N, l, K, L2, L3, J, I N V 

01000 

REAL  PI 

01100 

COMPLEX  X(0:  1023)/U,VI,T 

01200 

N  =  I  NT  (AEO'*fV  n«*l  .0  >  /  A  T  Or(7.0)  +  0.5) 

01300 

1F( ALOG(NUM*1.0)/ALOG(2.0)-N. GE. 0.001 ) 

GO 

TO 

3000 

01400 

J  =0 

01500 

PO  300  T=0,*!»iv 

OlbOO 

IF(I.GE.J)  GO  TO  100 

0-1 7  Ou 

T=v( J) 

OltiOO 

X(J)=X(I) 

01900 

v(n-r 

02000 

100 

K=(NUM+l)/2 

02100 

200 

TF(J.LT.K.PP.I.FG.NUM)  GO  TO  300 

02200 

J-J-K 

02300 

'(-Yll 

02400 

GO  TO  200 

02500 

300 

J  =.J  +  K 

02600 

F  1  =  3.1  41  *n->f 

02700 

OU  2000  1=1/ N 

02600 

I =2**1 

02900 

t.2 -L/2 

03000 

E3  -  L2-1 

03100 

”-(1.0/ 0.0) 

03200 

W  =  rMnLv(PCF(cI/I.2)/-SIN(ri/L2>) 

03300 

IF(INV.EQ.-I)  W=(1.0,0.0)/W 

0340O 

no  2000  J=0,L3 

03500 

no  1000  X=J,N"M,L 

03600 

T-X(K+L2)*U 

03700 

v(v+t  2 )  =  y ( ^ )-t 

03600 

100  0 

X(  f)  =  XCK)<-T 

03900 

2000 

U=U*R 

04000 

RETURN 

04100 

300C 

2 FITE  *OnO  3 

0  4*00 

4000 

FORMAT!'  THE  S'JMPER  OF  ELEMENTS  IN  THE 

INPUT' 

//, 

04200 

1  '  VECTCR  IS  NOT  2**N') 

04400 

RETURN 

04500 

END 
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00100 

o' 

ru  u;;k  suuvnu tines  in  in-:  i.ikrarv  run.  r-i;» 

00200 

C 

IN  A  nw  T»K  T  T.U1NJX  4006  S  FORAGE 

00300 

c 

SrU  PE  !•' V t-: ^ M ^  F  T"F  pnri'*  AM  Ar  cPLLOW^ 

00400 

c 

PTAKE  SYS:TEK6 

00500 

c 

00600 

c 

pEXErilTE  PI. GT.FOP,  PROGRAM  .  FOP,  SYS  :  TEKT.  KrL 

00700 

c 

00600 

c 

00900 

SHPPO»’TTMF  nLCT2  P ( X, Y, F ) 

01000 

c 

THIS  SUBROUTINE  PLOTS  T«E  VALUES  IN  THE  POSITIVE 

01100 

c 

REAL  2-n  A  A V  P(0:X,0:v).  T"E  PLOT  IS  MADF  UP 

01200 

c 

Y  +  i  H0RI70NTAL  LINES  EACH  CONTAINING  X+l  POINTS. 

01300 

c 

THE  MAXIMUM  HUMBER  OF  LINES  AND  POINTS  IS  512. 

01400 

c 

01500 

INTEGER  X,Y,  I,J,MINX,MINY,N,M,H 

01600 

HEAL  P(0:v,0:V),YMAX(0: 51 1),XINC,YINC,PMAV, RATIO 

01700 

WRITE  ( 5 , c0) 

01600 

50 

FORMAT ( *  TF  CnpY  IS  WANTED  TYPE  1  IF  NOT  TYPE  0' 

01500 

1  *  A^TEO  rl.0T  IS  FINISHED',/, 

02000 

2  '  WHpN  RPADY  SPACE  AND  RETURN') 

02100 

rFftn  ( 5,751  T 

02200 

75 

FORMAT!  Ii) 

02300 

Y  I!JC=51 2 . 0/Y 

02400 

R  AT  I 0  =  Ff  OAT( XVFL0AT(v) 

02500 

X1NC=512.0/V 

02t00 

FMAX=P(0,0) 

02700 

DO  100  1  =  0, X 

02800 

PO  100  -J-O/Y 

02900 

100 

IC(PMAX.LT.P(I/J))  PFAX=F(I,J) 

03000 

CO  200  T  =0,X 

03100 

200 

vf<Axm=o.p 

03*00 

CALL  INI TT (460) 

03300 

CALL  PWrNnQ(0.0»FLCAT(y.  )/0.Q,rpAX) 

03400 

DO  *00  •1  =  rt/v 

03500 

MI N  X=5t  2-XI NC* J 

03600 

MINY-YTNC*J 

03700 

CALL  SW IN DP (MI  NX, 511, MINY, 511  ) 

03800 

•'ALL  "nvpA(0. 0,0.0) 

03900 

DO  400  1=0, X 

04000 

IF  (P(T,J).LT.YMAX(I))  GO  TO  300 

04100 

CALI  r>FAWA(cLr’AT(n,P(I,J)) 

04200 

YMAX(T)=P(I, J) 

04300 

GH  TC  400 

04400 

300 

CALL  !-‘CVFA(FLOAT(I),P(I,J)) 

04500 

400 

r’nN'rIm,r 

04o00 

M=T NT (RATIO* FLOAT! J«-l ) )-I NT (RATI 0*FL0AT(J)) 

04700 

N=  v-|| 

04800 

no  500  1=0, N 

04900 

500 

y«  AY(  v_  T  )  =  YVAV(  )-P*'AX/Y 

05000 

M  =  M  -  1 

05100 

r0  600  I - 0 , F 

05200 

600 

YMAX( I )=P(  I  , J ♦  1  ) 

05300 

CALL  ',li''Ar,S(0,*'12') 

05400 

CALL  npvA6G( 512,0) 

05500 

CALL  ^RVIA^SI  1023,0  ) 

05600 

CALL  MOVAB3(0,767) 

05700 

CALL  AM"nn5 

05800 

RE An  (6,*)  H 

0550tj 

I  c  (  P .  F.Q  .  1  )  CALL  HrCOPV 

06000 

RETURN 

OolOu  t-  *' P 

06200  C 
ObJOO  C 

06  400  S!Hpi)»T!NF  "  L^Ti  n  f  •  ,  v*' !  N, \*>  A  X  ) 

OobOU  C  1-LiJTo  T’l  E  RE  A  L  ARRAY  P'OrNUM). 

06600  C  M'M  MAX  IS  1023. 

06700  C  YM1N  AND  YMA  X  ARK  THE  MINIMUM  AND  MAXIMUM 

06800  C  V  A  L"F.  9  nc  T"F  ApRAY  P.  T»'FV  ARP  POP  SpALTNr. 

36900  C  IF  T'.'U  PLOTS  A RF  TO  HAVE  THE  SAME  SCALE 

07000  C  THEM  THEY  M"ST  "AVF  THE  SAMF  VALUES  CF  vtfAX 

07100  C  AND  YMIN.  IF  T!)P  SUBROUTINE  IS  TO  DO  THE 

07200  C  SCALING  T”E"  pp.t  vvax  EpUAL  TO  VNIN. 

07300  C 

07400  INTEGER  I,NUM,8 

07S00  REAL  P ( 0 : NUM ) , YM I M, YM AX, XM I N, X MAX 

07600  /SITE  (C,',5M 

07700  2b  0  F  U  P  M  A  T  ( '  IF  COPY  IS  WANTED  TYPE  1  IF  NOT  TYPE  0',/, 

07800  1  '  A  c  T  E  0  PLCT  IS  PINISHEC',/, 

07S0O  2  '  WMF N  RF  A  p  Y  SPApE  AND  RETURN') 

08000  READ  ( b, 275)  I 

08100  27b  FORMAT (II) 

08200  IF(YMIN. NF.YMAX  )  CO  TO  350 

08300  Y V  I  R  =  ° ( 0  ) 

08400  YMAV-p(U) 

08500  DO  300  1=0, NUN 

08600  IF(pm.f  T.Y'-'IN)  YMTN=P(T) 

08700  300  IF(P(I).DT.YMAX)  YMAX-P(I) 

08800  Jbo  X»MN-:0.0 

08900  X!:AX=FLC  A?(NllM) 

09C00  r  ALT  IN  ITT  f A  pn  t 

09  100  CALL  ni.rNr'J(  XMIN/XMAX/YMIN/YMAX) 

09200  CALL  "C VE  A (0 . 0 , r ( 0  )  ) 

09  3  00  L'C  TOO  I  =  1  / *1 U M 

09  400  400  pALL  n®A<*ACFLrATf!)/p(T)) 

0 9 b 0 0  CALL  MUVAHSC  0/767) 

09t00  CALL  AN**GT'S 

09700  RFAP  (5,*)  H 

09800  P  FT” RN 

09900  END 


MISSION 

of 

Rome  Air  Development  Center 

RAVC  plan*  and  execute*  xeAeaxch,  development,  tut  and 
detected  acquisition  program*  in  *upport  o$  Command,  Control 
Communication*  and  Intelligence  (C3 1)  activitiu.  Technical 
and  engineering  *upport  within  area*  o&  technical  competence 
i*  provided  to  E SD  Program  OHicu  (PO*)  and  other  ESP 
element *.  The  principal  technical  minion  area*  axe 
communication*,  electromagnetic  guidance  and  control,  *ux- 
veillance  oh  ground  and  aero* pace  object*,  intelligence  data 
collection  and  handling,  information  *y*tem  technology, 
iono*phexic  propagation,  *olid  Atate  Adencu,  microwave 
pkyAic*  and  electronic  reliability,  maintainability  and 
compatibility. 


